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MULTIVARIATE ORTHOGONAL LAURENT POLYNOMIALS 
AND INTEGRABLE SYSTEMS 

GERARDO ARIZNABARRETA AND MANUEL MANAS 


Abstract. An ordering for Laurent polynomials in the algebraic torus (C*) D , inspired by the Cantero-Moral- 
Velazquez approach to orthogonal Laurent polynomials in the unit circle, leads to the construction of a moment 
matrix for a given Borel measure in the unit torus T°. The Gauss-Borel factorization of this moment matrix 
allows for the construction of multivariate biorthogonal Laurent polynomials in the unit torus which can be ex¬ 
pressed as last quasi-determinants of bordered truncations of the moment matrix. The associated second kind 
functions are expressed in terms of the Fourier series of the given measure. Persymmetries and partial persym- 
metries of the moment matrix are studied and Cauchy integral representations of the second kind functions are 
found as well as Plemej type formulae. Spectral matrices give string equations for the moment matrix which 
model the three-term relations as well as the Christoffel-Darboux formulae. 

Christoffel type perturbations of the measure given by the multiplication by Laurent polynomials are stud¬ 
ied. Sample matrices on poised sets of nodes, which belong to the algebraic hypersurface of the perturbing 
Laurent polynomial, are used for the finding of a Christoffel formula that expresses the perturbed orthogonal 
Laurent polynomials in terms of a last quasi-determinant of a bordered sample matrix constructed in terms 
of the original orthogonal Laurent polynomials. Poised sets exist only for nice Laurent polynomials which 
are analyzed from the perspective of Newton polytopes and tropical geometry. Then, an algebraic geometri¬ 
cal characterization of nice Laurent polynomial perturbation and poised sets is given; full column rankness of 
the corresponding multivariate Laurent-Vandermonde matrices and a product of different prime nice Laurent 
polynomials leads to such sets. Some examples are constructed in terms of perturbations of the Lebesgue-Haar 
measure. 

Discrete and continuous deformations of the measure lead to a Toda type integrable hierarchy, being the 
corresponding flows described through Lax and Zakharov-Shabat equations; bilinear equations and vertex op¬ 
erators are found. Varying size matrix nonlinear partial difference and differential equations of the 2D Toda 
lattice type are shown to be solved by matrix coefficients of the multivariate orthogonal polynomials. The dis¬ 
crete flows are connected with a Gauss-Borel factorization of the Jacobi type matrices and its quasi-determinants 
alow for expressions for the multivariate orthogonal polynomials in terms of shifted quasi-tau matrices, which 
generalize those that relate the Baker functions with ratios of Miwa shifted T-functions in the ID scenario. It is 
shown that the discrete and continuous flows are deeply connected and determine nonlinear partial difference- 
differential equations that involve only one site in the integrable lattice behaving as a Kadomstev-Petviashvili 
type system. 
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1. Introduction 

In this paper we give a theory for multivariate orthogonal Laurent polynomials in the unit torus (MVOLPUT). 
Our motivation comes from our previous work in orthogonal polynomials in the unit circle (OLPUC) 10, its 
matrix extension IfTTII and our recent developments on multivariate orthogonal polynomials in real spaces 
(MVOPR) fl3l ii4|. The natural framework for the study of the ring of Laurent polynomials in several vari¬ 
ables is the algebraic torus (C*) D , C := C \ (0). The orthogonal elements appear once a Borel measure d p 
with support in the D-dimensional unit torus T°, T := (z G C : |z| = 1}, is considered. Inspired by the 
Cantero-Moral-Velaquez (CMV) approach [23' [94), se also EHII, we propose an order for the multivari¬ 
ate Laurent polynomials in the algebraic torus (C*) D . For the MVOLPUT we find in this paper not only 
three-term relations and Christoffel-Darboux formulae, but also construct associated second kind functions 
and partial second kind functions given in terms of the Fourier series of the measure and Cauchy integral 
transforms, all connected through Plemej relations. Remarkably, we are able to find the extension of the 
Christoffel formula for Laurent polynomial perturbations of the measure in terms of a last quasidetermi¬ 
nant of a bordered sample matrix constructed from the original MVOLPUT evaluated at a poised set in the 
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algebraic hypersurface of the zeroes of the Laurent polynomial generating the perturbation. As the chosen 
order does not generate any gradation, for the finding of poised sets we need of what we have called nice 
Laurent polynomials. The discussion requires of mathematical ideas connected with tropical geometry as 
Newton polytopes and also of elementary algebraic geometry. Finally, we are able to construct Toda inte- 
grable systems of discrete and continuos type and also find how the Kadomtsev-Petviashvilli flows live in 
each site of the Toda lattice. 

For a better understanding of the results achieved and its context let us perform some preliminary com¬ 
ments on the state of the art of the subjects involved. 

1.1. Historical background. 


1.1.1. Orthogo7ial polynomials in the unit circle. The unit circle T is the border of the unit disk D := {z 6 C : 
\z\ < 1}. A complex Borel measure p supported in T is said to be positive definite if it maps measurable 
sets into non-negative numbers, that in the absolutely continuous situation (with respect to the Lebesgue- 
Haar measure has the form w(0) d 0, 0 6 [0, 2 tt)]. For the positive definite situation the orthogonal 
polynomials in the unit circle (OPUC) or Szego polynomials are defined as the monic polynomials P n 
of degree n that satisfy the orthogonality relations, J T P n (z)z k d pfz) = 0 , for k = 0 , 1,..., n — 1, 
Orthogonal polynomials on the real line (OPRL) with support on [—1,1] are connected with OPUC 
In the real case the three-term relations provide a tridiagonal matrix, the so called Jacobi operator, while 
in the unit circle support case the problem leads to a Hessenberg matrix. OPUC's recursion relation is 
given in terms of reverse Szego polynomials P{(z) := z l Pdz ' ) and reflection or Verblunsky coefficients 
(Xi := Pi(0). The recursion relations for the Szego polynomials can be written as 


' Z OCI 

, Z Oil 1 


) (p*_] 


Szego's theorem implies for a nontrivial probability measure d p on T with Verblunsky coefficients {a n }£V 0 
that the corresponding Szego's polynomials are dense in L 2 (T, p) if and only if — |a n )| 2 ) = 0- For an 

absolutely continuous probability measure Kolmogorov's density theorem ensures that density in L 2 (T, p) 
of the OPUC holds if and only if the so called Szego's condition J T log(w(0) d 0 = —oo is fulfilled, It88l . We 
refer the reader to Barry Simon's books |85| and If86l for a very detailed studied of OPUC. 


1.1.2. Orthogonal Laurent polynomials and CMV. In the context of the strong Stieltjes moment problem the 
papers |[57i |58 l could be considered as the seed for the study of orthogonal Laurent polynomials on the 
real line (OLPRL). Whenever we have solution of the moment problem we have Laurent polynomials 
{Qnirvdo that fulfill the orthogonality conditions J R x~~ n+ ^Q n (x) d pfx) = 0 for j = 0,... ,n — 1. For the 
early development of the theory see [24;, 33j [56| and Il75| . The circle T context was first considered in M90l 
, see also Il28l |2 3i [26l 271 where matters like recursion relations, Favard's theorem, quadrature problems, 
and Christoffel-Darboux formulae were treated. The CMV If23l representation is a hallmark in the study of 
certain aspects of Szego polynomials. Indeed, despite the set of OLPUC being dense in L 2 (T, p) in general 
this is not true for the OPUC, Il22l and [28l . Now, the recursion relations for ordinary Szego polynomials 
and its reverse polynomials is replaced by a five-term relation similar to the OPRL situation. Alternative 
or generic orders in the base used to span the space of OLPUC can be found in (27]. Other papers have 
reviewed and broadened the study of CMV matrices, see for example (871I6IH . As was pointed out in l87l 
the discovery of the advantages of the CMV ordering goes back to previous work If94l . 

1.1.3. Midtivariate orthogonal polynomials. The monographs |37| and 11011 are highly recommend references 
for understanding the state of the art regarding multivariate orthogonal polynomials. The recurrence re¬ 
lation for orthogonal polynomials in several variables was studied by Xu in l96l . while in ||97l he linked 
multivariate orthogonal polynomials with a commutative family of self-adjoint operators and the spec¬ 
tral theorem was used to show the existence of a three-term relation for the orthogonal polynomials. He 
discusses in |98l how the three-term relation leads to the construction of multivariate orthogonal polyno¬ 
mials and cubature formulae. The analysis of orthogonal polynomials and cubature formulae on the unit 
ball, the standard simplex, and the unit sphere 11001 lead to conclude the strong connection of orthogonal 
structures and cubature formulae for these three regions. The paper (99| presents a systematic study of 
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the common zeros of polynomials in several variables which are related to higher dimensional quadra¬ 
ture. Karlin and McGregor |[60l and Milch Ii7ll discussed interesting examples of multivariate Hahn and 
Krawtchouk polynomials related to growth birth and death processes. There have been substantial de¬ 
velopments since 1975, for instance, the spectral properties of these multivariate Hahn and Krawtchouk 
polynomials have been studied in H47l . In l46l a two-variable positive extension problem for trigonomet¬ 
ric polynomials was discussed —here the extension is required to be the reciprocal of the absolute value 
squared of a stable polynomial. This could be understood as an autoregressive filter design problem for 
bivariate stochastic processes. The authors show that the existence of a solution is equivalent to solving a 
finite positive definite matrix completion problem where the completion is required to satisfy an additional 
low rank condition. As a corollary of the main result a necessary and sufficient condition for the existence 
of a spectral Pejer-Riesz factorization of a strictly positive two-variable trigonometric polynomial is given 
in terms of the Fourier coefficients of its reciprocal. A spectral matching result is obtained, as well as in¬ 
verse formulas for doubly-indexed Toeplitz matrices. Geronimo and Woerdeman used tools including a 
specific two-variable Kronecker theorem based on certain elements from algebraic geometry, as well as a 
two-variable Christoffel-Darboux like formula. In Il62l reproducing kernels are used to give and alternative 
proofs for some of the mentioned results on orthogonal polynomials on the two dimensional torus (and 
related subjects). Regarding extension problems on the torus and moment matrices see Ifl5l . In I134H a for¬ 
mula describing the asymptotics of a class of multivariate orthogonal polynomials with hyperoctahedral 
symmetry as the degree tends to infinity is given. The polynomials under consideration are characterized 
by a factorized weight function satisfying certain analyticity assumptions. 

1.1.4. Darboux transformations. In the context of the Sturm-Liouville theory, Gaston Darboux Il29l intro¬ 
duced these transformations for the first time. Much later |[70| this transformation was named after Dar¬ 
boux. It has been extensively developed for orthogonal polynomials 11501 [95 , 201 2B E9|. In Geometry, the 
theory of transformations of surfaces preserving some given properties conforms a classical subject, in the 
list of such transformations given in the classical treatise by Einsehart [[38] we find the Levy (named af¬ 
ter Lucien Levy) transformation l64l . which later on was named as elementary Darboux transformation 
and known in the orthogonal polynomials context as Christoffel transformation Il95l [841. The adjoint ele¬ 
mentary Darboux or adjoint Levy transformation is also relevant Il70l [35 1 and is referred some times as a 
Geronimus transformation Il95l . Linally the rational Uvarov transformation corresponds to the fundamen¬ 
tal transformation introduced in a geometrical context by Hans Jonas 155] , see also the discussion in |38| . 
Lor further information see [78J, 51|. The iteration formula of Christoffel transformations is due to Elwin 
Bruno Christoffel B25l . This fact was rediscovered much latter in the Toda context, see for example the 
formula (5.1.11) in fl7dl for W+(N). 

1.1.5. Integrable systems. The papers [81.. 821 and I1301I3T1I32| settled the Lie group theoretical description of 
integrable hierarchies. See |[74l for a discussion of the role of the factorization problems, dressing proce¬ 
dure, and linear systems as the keys for integrability. In this dressing setting the multicomponent integrable 
hierarchies of Toda type were analyzed in [91,192:. [93]. Lor further developments see lfl8llT9ll and [59.681 for 
the multi-component KP hierarchy and |67| for the multi-component Toda lattice hierarchy. Mark Adler 
and Pierre van Moerbeke showed how the Gauss-Borel factorization problem appears in the theory of the 
2D Toda hierarchy and what they called the discrete KP hierarchy |lj !2K3 SISJ [6]. These papers clearly 
established -from a group-theoretical setup- why standard orthogonality of polynomials and integrability 
of nonlinear equations of Toda type where so close. In fact, the Gauss-Borel factorization of the moment 
matrix may be understood as the Gauss-Borel factorization of the initial condition for the integrable hier¬ 
archy. 

1.1.6. Our previous work. In the Madrid group, based on the Gauss-Borel factorization, we have been 
searching further the deep links between the Theory of Orthogonal Polynomials and the Theory of In¬ 
tegrable Systems. In j9] we studied the generalized orthogonal polynomials HU and its matrix extensions 
from the Gauss-Borel view point. In IfTOl we gave a complete study in terms of this factorization for mul¬ 
tiple orthogonal polynomials of mixed type and characterized the integrable systems associated to them. 
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Then, we studied Laurent orthogonal polynomials in the unit circle trough the CMV approach in J7J and 
found in |j8]| the Christoffel-Darboux formula for generalized orthogonal matrix polynomials. These meth¬ 
ods where further extended, for example we gave an alternative Christoffel-Darboux formula for mixed 
multiple orthogonal polynomials Ifl2l or developed the corresponding theory of matrix Laurent orthogonal 
polynomials in the unit circle and its associated Toda type hierarchy IITT1 . 

1.2. Layout of the paper. In §2, for a Borel measure in the unit torus this monomial ordering induces an 
arranging of the moments in a moment matrix, built up of growing size rectangular blocks, as it happens for 
MVOPR ii, this moment matrix has a rich structure having therefore several important properties. In 
the quasi-definite case the Gauss-Borel factorization of this matrix gives multivariate biorthogonal Laurent 
polynomials in the unit torus and to MVOLPUT for positive measures. We define second kind functions 
that are shown to be a product of the MVOLPUT and the Fourier series of the measure. The moment 
matrix has persymmetries and partial persymmetries which happen to be useful for the finding of Plemej 
type formulae for the orthogonal Laurent polynomials. Then, another symmetry of the moment matrix, 
which we call string equation, leads to three-term relations and Christoffel-Darboux formulae. 

We study in §3 Laurent polynomial perturbations of the measure, which could be considered as an 
extension of the Christoffel transformation to this framework. Interpolation theory and the construction 
of appropriate sample matrices based on poised sets allows for the construction of the Christoffel formula 
in this generalized scenario. To construct poised sets we need of what we call nice Laurent polynomials, 
which behave nicely with the extended CMV ordering. We study nice Laurent polynomials with the aid 
of Newton polytopes, and some aspects of tropical geometry. Full column rank Laurent-Vandermonde 
matrices and algebraic geometry lead to poised sets in this context. These Darboux transformations are 
studied for the Lebesgue-Haar measure. Perturbed Christoffel-Darboux kernels and connection formulae 
for the perturbed and non perturbed kernels are given 

Finally, in §4 we introduce integrable systems of Toda and Kadomtsev-Petviasvilii (KP )type. In par¬ 
ticular we show how discrete Toda type flows appear and give Miwa type expressions for the orthogo¬ 
nal Laurent polynomials. A Toda type integrable hierarchy is found, and Baker functions. Lax matrices 
and Zakharov-Sahabat equations are discussed. We also find bilinear equations, Miwa shifts and corre¬ 
sponding vertex operators. Finally, the use of asymptotic modules allows for the construction of KP onsite 
flows; i.e. to integrable nonlinear equations constructed in terms of the coefficient associated to a given 
MVOLPUT, say cf>[kl with a fixed k. 


2. Multivariate biorthogonal Laurent polynomials in the unit torus 

The algebraic torus (C*) D , C* := C \ {0}, is an Abelian group under component-wise multiplication and 
has as its coordinate ring the ring of Laurent polynomials C[z ±1 1 = C [zf 1 ..., Zq 1 ], where we require of D 
independent complex variables z = (zi, Z 2 ,..., zd) T £ (C*) D . Therefore, the algebraic torus is the natural 
framework when considering the ring of multivariate Laurent polynomials. 

D 

Given the unit torus T D = X T, the Cartesian product of D copies of the unit circle, and a a Borel measure 

i=l 

d p £ B (O) on some complex domain O such that its support belongs to the unit torus, supp p C T°, we 
will study the corresponding orthogonal Laurent polynomials and their properties. Remember that the 
D-dimensional unit polydisk D° has the unit torus T D as its distinguished or Shilov border. For z £ T d we 
have the parametrization 

z(0) = (e i0 \...,e i0D ) T 


with 0 := (0i,..., 0 d ) t £ [0,27t) D . A complex Borel measure p supported in T D is said to be real, re¬ 
spectively positive definite, if it maps measurable sets into the real numbers, respectively into non neg¬ 
ative numbers. When the measure is absolutely continuous with respect to the Lebesgue-Haar measure 
d0 = d0i • • • d 0 d it has the form dp(0) = w(0) d 0 and the weight w(0) must be a positive function when 
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p. is positive definite. The inner product of two complex valued functions f (z) and g(z) is defined by 


(f, g) := o f(z(0))g(z(0)) d p(0). 
Jt d 


2.1. Laurent monomials. Ordering Laurent polynomials. The ring of complex Laurent polynomials 

is the localization of the polynomial ring C[zi,..., Zq] by adding the formal inverses of z \,..., Zd- The 
units; i.e v invertible elements in the ring, are the Laurent monomials az“, a 6 C*. Following Il42l we say 

Definition 2.1. A Laurent monomial is a well-defined function between algebraic tori z“ : (C*) D —* C*. 

A Laurent monomial can be viewed as a character of the algebraic torus and Laurent polynomials can 
be thought as finite linear combination of Laurent monomials. Any usual polynomial in the ring C[z] can 
be considered as a Laurent polynomial. 

Definition 2.2. (1) Given a multi-index a = (ai,..., cxd) t £ Z D we write z a = z“' • • • . 

(2) For a finite subset A C Z D we define C- 4 as the set of Laurent polynomials with monomials from A; L £ C A 
if and only if L = Y. L a z“ with L a f 0 for a £ A. 

(3) When L £ C- 4 we say that A is the support of L and its convex hull is knoivn as its Newton polytope NP(L) 

NP(L) := Conv(A). 

(4) A supporting hyperplane H of the Nezvton polytope NP(L) is an affine hyperplane such that the face F = 
H n NP(L) A 0 an & NP(L) is fully contained in one of the two halfspaces defined by H. A face is called facet 
if it has codimension 1. Given w £ M D we construct the associated face 

F W (NP(L)) := {u £ NP(L) : (u-v) • w < 0 Vv £ NP(L)}. 

Given a face F o/NP(L) its normal cone is A/"f(NP(L)) := {w £ M D : F = F W (NP(L))}, the normal fan 
A/"(NP(L)) := |A/r(NP(L)) : F is a face o/NP(L)| is the collection of all normal cones. 

(5) For any multi-index a £ Z D we define its longitude |a| := ]Tq=i |a a |. 

(6) Given a Laurent polynomial L £ C A its longitude is 

£(L) := max |a|. 
aeA 

Observe that for the ring of polynomials C[z] this longitude is the total degree f37, 131, and given two 
polynomials Pi, P 2 £ C[z] we have deg(Pi P 2 ) = deg Pi +deg P 2 Q However, for the Laurent polynomial ring 
we have f(Li L 2 ) < i (Li) + T(L 2 ) j^] Notice that this order is not a monomial order; i.e., given two monomials 
Li < L 2 then the multiplication by any other monomial L does not respect in general this order and it could 
happen that LLi > LL 2 . Despite the previous observation we should stress that the corresponding lattice is 
a graded lattice. Our proposal is motivated by the D = 1 case where one has the CMV ordering |23j[94j, for 
the Laurent polynomial ring Clz^ 1 ], in where powers z n and z ~' n go together 123[|Zt- Despite that there is Z- 
grading, deg z^ = ±1, we reckon that it is not as useful for the analysis of orthogonal Laurent polynomials 
as the CMV ordering is. This will be relevant later on when we discuss Darboux transformations. 

Definition 2.3. Given a nonnegative integer k £ Z + zve introduce 

[k] :={«£ Z D : |a| = k}. 


lr rhe ring of polynomials C[zi,..., zd] is a Z+-graded ring C[zi, ..., zd] = © nez+ C n [zi,..., zd] with C n [zi,..., zd] the 
degree n homogenous polynomials and C n [zi,... ,zo]C m [zi,... ,Zd] C C n +m[zi, • •. ,zd]- 

2 This follows from |n + m| < |n| + |m|, that holds for all n, m G Z, inequality that is saturated |n + m| = |n| + |m| for either 
couples of positive integers or couples of negative integers. 
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Proposition 2.1. The number of multi-indices in TP of longitude k > 0 is 


( 2 . 1 ) 


min(k,D) 

I Ml = X 2* 

j=1 


D\ /k- 1 


j-i r 


and | [0] | = 1. 

Proof. Recall that there are (?) strings containing j ones and (D — j) zeros, and observe that, when counting 
compositions, the number of ways to write k = ai + ag + • • • + ay where every eg is a positive integer, 
is given by (p,)• Thus, (?) (Ip,) is the number of partitions k = ai + a .2 + • • • + Qd, with non-negative 
integers at E Z+ for i E (1,..., D}; i.e., of weak compositions, with (D—j) factors equal to zero and j positive 
integers. We now drop the non negative condition and allow for arbitrary integers, k = I ai | +1 a .2 1 + • • • +1 ao I 
as we have j non-vanishing components, we should multiply by 2, as we must take into account a reversal 
of sign, that gives the same longitude, and T ] (?) (p|) should be the number of multi-indices a E Z D of 
longitude k having j of its D components different from zero. Finally, summing up in j we get |[k]|. □ 

Definition 2.4 (The longilex order in Z D ). |^Fz>sf, we order according the longitude of the multi-indices; i.e., 
a < a' whenever |ct| < |a'| Second, we use the lexicographic order for the set M of multi-indices of same longitude 
and write 



(X 


0c) 



< a 


(k) 

a+1' 


Proposition 2.2. (1) The convex hull Conv(M) is a regular hyper-octahedron with vertices given by 


{V| k ]:= (±k, 0,..., 0) 1 , Vp) := (0, ±k, 0,..., 0) 1 ,..., V£ J ± := (0,..., 0, ±k) 1 } C 


iT 


/(k) 


(k) 


(2) The Newton polytope of a Laurent polynomial belongs to the regular hyper-octahedron NP(L) C Conv([£(L)]), 
and at least a face of the Newton polytope NP(L) has a nontrivial intersection with a face o/Conv([£(L)]). 


Definition 2.5. We introduce the semi-infinite vector x constructed by using the longilex order of Laurent 
monomials 


( X[o]\ 
X[i] 


X := 


X[k] 

V = 


where 


X[k] 


( z“' \ 
z a 2 


\ Z a IMI 


In particular, the two first are easy to write 


[k] = { 


a 


(k) „(k) 

1 






(k) 

Ilk] I 


}• 


( 2 . 2 ) 


X[0] = 1 > 



X[i] 


.-1 


zd 




Z2 

\Zl / 


3 This not the deglex ordering introduced in the SAGE package for the Laurent ring polynomial. 
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and for D = 2,3 we have 


/ z 1 2 \ 


X[2] = 


Z 1 z 2 

V Z2 

—2 


Zi Z2 
ZtZ 2 


X[2] = 


Z 1 lz 2 1 

V 1 

1, 


Z 1 z 3 


•‘-I Z3 

zf 1 z 2 

z 2 2 

Z 2 1z 3 

£ 

z 3 

z 2%' 

Z2Z3 


ziz 2 

Z 1 Z 3 1 

ZlZ3 

ZiZ 2 


V 




respectively. 

In this paper we will consider semi-infinite matrices A with a block or partitioned structure induced by 
the longilex order of Definition |2 .4 1 

\ 


A = 


( A [0],[0] 

A [o],m • • A 


^ A «' k) ,a< l) ’ 

A (k) 

a. 

A m,[o] 

A m,m • • • 

, A [k],[l] = 



V i 

: J 


t A a (k) a (l) • 

• A a ,k) 
a IMI 


G 


^<|[k]|x|[l]| 


We use the notation G C lk |x l[i] f or tp e rectangular zero matrix, 0- k ; G C A for the zero vector, and 
I[ k ] G C ' Ax A £ or £p e identity matrix. For the sake of simplicity and if there is no confusion we prefer to 
write 0 or I for the zero or identity matrices, and we implicitly assume that the sizes of these matrices are 
the ones indicated by its position in the partitioned matrix. 


Proposition 2.3. The linear space C^[z ±] ] of Laurent polynomials of longitude k has dimension := dimC^lz 
given by 


±11 


N k = 


x=o 


2.2. The moment matrix. 


Definition 2.6. Given a Borel measure p with support in the unit torus T° the corresponding moment matrix 
is given by 


(2.3) 


G h :=d> 


T D 


x(z(0))dji(0)x(z(-0)) 


For the sake of simplicity when not needed we omit the subscript p and write G instead of G^. We write the moment 
matrix in block form 


/Gn 


G = 


Wo],m 
; m,[o] G [1]}[1] 
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with each entry being a rectangular matrix with complex coefficients 


(2-4) G [lc])[l] :=o X[k](z(9))d^0)X[i](z(-0)) 

Jt d 


/ G (ic) (i) 


< (k) a ,l) 

*1 >“ m 


\ 


V G a ,k) a ,l) G a (k) a (l) / 


G C IM|x|[I]| 


J (k) (l) • — 

ot, ,« 


4> 


T D 


■i (k) (l), „ 

e l( «^ ~°b ) ' 0 d|T(0)GC. 


Truncated moment matrices are given by 


G [k] := 


( G L0],[0] 

\G[k—!],[( 


G[o],[ic—i] 
G[lc—1],[k—1] 


and for l ^ k we will also use the following bordered truncated moment matrix 


r [k] . 
b l — 


/ G [0 ],[o] ••• G [ 0 ],[k-i] \ 

G [lc—2],[0] ••• G[k_2],[k-1] 

V G[1],[0] ••• G^j^.u J 


where we have replaced the last row of blocks, (G^jo] • • • G [ic_i],[ic_i]) ^ of the truncated moment matrix G^ by 
the row of blocks (G^jq] • • • G^^k-i])/ we also need a similar matrix but replacing rows by columns 


/ G [0 ],[ 0 ] 

' G [0]> [k—2] 

G [o],[i] \ 

V G[ic~i],[o] • 

• G [k —i ] 5 [ k — 2 ] 

G[k—1],[L] / 


Let us extend to this scenario the concept of quasi-definite and positive definite 
Definition 2.7. The measure and its moment matrix are quasi-definite if all its principal block minors are not singular 

detG [k] /0> k e {0,1,...}. 


When all the minors are positive 

detG [k] > 0, 


Ice {0,1,...}. 


we say that it is a definite positive moment matrix 

Notice, that instead of positive definite we could request the moment matrix to be definite; i.e., all block 
principal minors are positive or all are negative. 

We are now ready to discuss some aspects regarding the Gauss-Borel factorization of this moment ma¬ 
trix. 


Proposition 2.4. (1) A quasi-definite moment matrix G admits the following block Gauss-Borel factorization 

(2.5) G = S~ 1 H(§^ 1 ) t , 
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with 


( %0]| 

0 

0 

...\ 


%i]| 

0 


(S _1 )[2],[0] 

( S_1 )[2],[1] 

I |[2]i 


V : 




/H [0] 0 

0 • • A 



0 H m 

0 



0 0 

H[2] 



V : : 




( %]| 

0 

0 

..A 


%1]l 

0 



(S-’Jkui] 

h2]i 


V : 





Moreover, the quasi-tau matrices H[ k j are non singular, k € {0, 1 ,...}, and 

i-i 

det G ^ det H[ k ]. 

k=0 

(2) The measure is real if and only if its moment matrix is Hermitian G = G f In this case 

§ = S, H f = H. 


(3) The measure is definite positive if and only if the moment matrix is definite positive. In this case, the quasi-tau 
matrices H[ k ] are definite positive matrices for k e {0,1,...}. 


Proof. See Appendix C.l 


□ 


Quasi-determinants are extensions of determinants —more appropriately of quotient of determinants—> 
and fulfill the heredity principle, quasi-determinants of quasi-determinants are quasi-determinants. The 
Gel'fand school has given a very complete study on the subject, see Il42ll43ll44ll45l . However, in this paper 
we require of the generalization given by Olver in |[76l . A last quasi-determinant version of the above result 
can be given 


Proposition 2.5. If the last quasi-determinants of the truncated moment matrices are invertible 

det©*(G [k] ) /0, k = 1,2,... 

the Gauss-Borel factorization ( |2.5| can be performed where 

H [k] =0*(G [k+1] ), (S- 1 ) [k])[l ] = 0*(G^' +1] )0*(G [l+1] )- 1 , (S" 1 )^ = (©^G^V®*^'”))*. 

Proof. It is just a consequence of Theorem 3 of Il76ll . 

Definition 2.8. The matrices H[ k ] are called quasi-tau matrices. We introduce the first subdiagonal matrices 


P[k] : = 

S[1c],[Tc—1], 

$[k] ■= S[k],[k-1], 




k > 1, 


which take values in the linear space of rectangular matrices C^ x ^ k G 

and also define 




( 0 

o 

o 

o 


/ o 

0 

0 

0 • 

•\ 


P[i] 

0 0 0 ••• 


^[i] 

0 

0 

0 • 


(3 = 

0 

|3[2] 0 0 

, 0 = 

0 

^[2] 

0 

0 • 


0 

0 |3 [3] 0 ••• 

0 

0 

£[3] 

0 • 



V : 

; •. ••• ■■■) 


V ! 
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An immediate consequence of Proposition |2.5|is 







□ 
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Proposition 2.6. The first subdiagonal matrices have the following quasi-determinantal expressions 

Pm ^-0*(gJ c1 )0*(g m )- 1 , = -(©*(6f)e*(G M r 1 ) t . 

2.3. Orthogonal Laurent polynomials in the unit torus. With the aid of the Gauss-Borel factorization we 
introduce 


Definition 2.9. We define the multivariate Laurent polynomials in D complex variables z \,..., zd 


( 2 . 6 ) 



( ^ 


k 

(^ \ 

O := Sx = 

4*11] 

> 4>[k]( z ) 

= y S[ k]| [ l ]X[i]( z ) = 



l ! ) 


1=0 



/$ [0 A 


k 


/N /\ 

O := Sx = 

$m 

> $[k] ( z ) 

= y ^[ki,[i]X[i]( z ) = 



(i J 


1=0 



1=0 j=l ’ 


(X) 


-.LX S a ! k V l,Z ' 

! 1=0 j=l ’ 




Observe that <J)[ k ](z) = X[k]( z ) + P:k]X[k-i] ( z ) + • • • is a vector constructed with the multivariate Laurent 
polynomials (z) of longitude k, each of which has only one monomial of longitude k; i. e., we can write 
4) ai (z) = z at + Q a t (z), with £(Q a J < k; the same holds for $ M . 

Proposition 2.7. The two sets of Laurent polynomials {<t>[k]}k°=o am ^ {$[k]}uLo f orm a biorthogonal system; i.e., 


(2.7) 


o 4^[ k] (z(0)) d p(0) ($ tl] (z(0))) f =6 k)l H [lc] . 
Jt d 


Proof. It is straightforward from the Gauss-Borel factorization. 


□ 


This biorthogonality leads to the following orthogonality relations 

o <t ) [ k] (z(0))d[x(0)(x [l ](z(0))) t =o X[i]( z (0))dp(0)($ [k] (z(0))) t = 0, l = 0,..., k — 1 , 

Jt d Jt d 

o 4> [k ](z(0)) d (j.(0)(x[ k ](z(0))) f = o X[k]( z (0))dp(0)($ [k ](z(0))) t = H [k] , 


JT d 


Jt d 


Proposition 2.8. When d p(0) is positive definite the set {cj)[ k ](z)}^ 0 is an orthogonal set of Laurent polynomials 
in the unit torus T D ; i.e., 


( 2 . 8 ) 


o 4^ [k ](z(0)) d |x(0)(4) [l] (z(0))) t = 5 k)l H [k] . 
Jt d 


For k = l we get 


o ||4 ) [k] (z(0))|| 2 d p(0) = H [k ], 
Jt d 


and therefore the quasi-tau matrices H[ k ] may be viewed as the squared norm matrices of the MVOLPUT. 
Notice that we are talking about biorthogonality or orthogonality in a block form. For example, for real 
measures what we have is that the blocks of different longitude are orthogonal, which is fine, however the 
Laurent polynomials of the same length are not truly orthogonal but satisfy 


cb <|> ( ki(z(0))dp(0)* (io(z(0)) = 6 k) iH (k) (k) . 

J T D “i “j “i 

In the positive definite case the quasi-tau matrix H[ k ] is a positive definite Hermitian matrix and we could 
find and orthogonal transformation giving standard orthogonal Laurent polynomials. But if we do so we 
spoil the symmetry that leads to fundamental properties, as we will see later, of these polynomials. 
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Proposition 2.9. The following quasi-determinantal formuleehold true 


/G[, 


4*[k] (z) — ©* 


••• G[ 0 ],[ic— 1 ] X[o](z) N 


($[ k ](z)) f = 0* 


\G[k],[ 0 ] ••• G [k],[k- 1 ] X[k](z), 


/ G[, 


Grv- 


[k—!],[( 


V(X[0]( z )) 


G[0],[k] \ 


G[k—1],[k] 
(X[k](z))V 


2.4. Holomorphic extensions of the Fourier series of the measure and second kind functions. First we 
recall some basic facts regarding the analysis in several complex variables, see flT6ll53l 83 :63. 791 for more 
information. Given the vector with positive components r = (ri,..., ro) T G M+, the polydisk 

D D (r) = {z = (z 1 ,...,z D ) T : |zt| < n,i G {1, ..., D}} C C D 

centered at the origin of polyradius r has as its distinguished boundary, also known as Shilov border, the 
D-dimensional torus 


T D (r) = {zGC D :|z i |=r i ,iG{l,...,D}}. 

For any two polyradii r and R the associated polyannulus centered at the origin is 

A D (r, R) := {z G C D : n < z t < Rpi G {1,..., D}}. 


Let us recall that a set A C C D is a complete Reinhardt domain if the unit polydisk O d acts on it by com¬ 
ponentwise multiplication. The polydisk of convergence of a power series is such that any other polydisk 
D d (r 4 5 6 7 ) with Tj < r- for some j G {1,..., D} contains points in where the power series diverge. 

We remind the reader that any set A C C D is a Reinhardt domairj^] if the unit torus T D acts on it (for 
every c G A and e 10 G T D we have that (e 101 Ci,..., e 10D Cq) t G A). The domain of convergence T> j_ C 
C D of a Laurent series L(z) = LaGa z“ is a Reinhardt domain. Recall that for all polyradii r and R the 
annulus A D (r, R) is a Reinhardt domain and that any Reinhardt domain is the union of polyannuli. The 
Laurent series is locally normally summable in its domain of convergence and therefore locally absolutely 
uniformly summablerl The function L(z) is holomorphic (holomorphic in each variable Z{, i G (1,..., D}) 
in T>i, which is its domain of holomorphy. Conversely, given a holomorphic function L(z) in AP q (r, R) (a 
polyannullus centered at Zq G C d ), and a polyradius p such that p < p t < R L , i G { 1,..., D}, then 


L(z) = Y_ l- a (z-z 0 ) a , 

aez D 


L„ = 


(27ti 


|D 


L(z) 


T D (z 0 ,p) [z-ZqY 


dzi . . . d Zq , 


where T D (zo, p) is the distinguished border of the polycircle centered at Zo with polyradius p. 

We consider the multivariate Fourier series (1(0) of the measure p given in terms of its moments or 
Fourier coefficients 


(2tt) 


D 


T d 


e^ ia ' e dp(0), 


cx GZ D , and (1(0) := ^ c a e 1<x ' e . 

aez D 


As c a (p) = C—cx)p), for real measures we deduce that = c a and, consequently, p(0) = (t(—0). Let 

D(T d ) ={f G C°°(M d ) :f(x + 27tcx) = f(x), Vx G M D , Vex G Z D } 


be the linear space of test functions. Then, the Fourier series always converges in D 7 (T d ), the space of 
distributions on the D-dimensional unit torus Il80l . so that J TD (t(0)f(0) d 0 = J TD f (0) d p(0), Vf G D(T d ). 
For an absolutely continuous measure d p(0) = w(0) d 0 we can write d p(0) = (1(0) d 0. If we assume 
that p is a Radon measure, i.e. locally finit^] and inner regular]^] then the associated linear functional 
f ha J tD f(0) d p(0) is continuous and therefore is a distribution. 


4 Sometimes named D-circled domain or circled domain. 

5 A Laurent series 2I aeZ D L a z“ is locally normally summable if for any compact set K C I>l there exists C > 0 and 0 € (0,1) 
such that |L a z“| < C0 |oi| for z G K and a £ Z D . 

6 V0 G T D exists a neighbourhood U such that p(U) < 00 . 

7 p(S) = sup KcB p(K) where the K are compact subsets of the Borel set B. 
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Given the Fourier coefficients or moments c a of the measure d p, generating its Fourier series (1(0), we 
have the corresponding Laurent series 

P-(z) := c “ z “ 

aeZ D 

that converges in a Reinhardt domain V^. This domain, which is a union of polyannulli, is the of domain 
of holomorphy of jl(z), and the series is locally absolutely uniformly summable there. Notice that the 
Reinhardt domain of the Laurent series of the measure belongs entirely to the algebraic torus V^ C (C*) D . 
In terms of the Fourier coefficients the moment matrix \2A) can be expressed as 


Motivated by Definition 2.6 we now consider 

Definition 2.10. The second kind functions are defined by 

C = (S- 1 )t x (z), (?=(g- 1 )t x (z). 

We now show how these functions are linked with the MVOLPUT 

Proposition 2.10. Let V ^ be the Reinhardt domain of convergence of the Laurent series jl(z). Then, the second kind 
functions can be written in terms of p.(z) and the MVOLPUT asfollozvs 


€ V 


in 


£(z) =(27t) D H _1 <h(z)(l(z), zeP p 


C(z) = (2 7 r) D (Ht)- i a(z)fi(z- 1 ), 

Proof. See Appendix C.2| □ 

2.5. Persymmetries. The moment matrix is a structured matrix. Indeed, by construction according to the 
longilex order, it is subject to important symmetries. Let us discuss now the persymmetry that it fulfills. 
We notice that if a G [k] then —a E [k]; moreover, using the longilex order of Definition |2.4| we find 


(k) 

‘iMI+l-i 


„(k) 


To model this fact we introduce 


Definition 2.11. Given any non negative integer m E Z + we consider the exchange matrix 

/0 0 
0 0 
0 0 

£m 

0 1 
V 0 

Then, zve introduce the follozving semi-infinite block diagonal reversal matrix 

n : = diag(%|,%]|,...). 

The exchange matrix is also known as reversal matrix, backward identity matrix, or standard involutory 
permutation matrix, see [ 52 . (491. 

Proposition 2.11. The follozving properties hold true 

t) =n _1 =r) f = r| T . 

The vector X fulfills 

hx(z) = x( z 1 )> (x(z)) f n = (x(z _1 )) t , 

zvhere z _1 := [zf ',..., zff ) T ; in components, the previous relation reads 

(k) (fc) 

(%]|X[k](z))i = Z^MI+i-i = Z~N . 


0 0 1 \ 
0 1 0 
1 0 0 


0 0 0 
0 0 0 / 




1, j =m — i + 1, 
0, j f m — i + 1. 
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In terms of the reversal matrix £ m a matrix M e C mxm is said to be persym metric, see Ii49i j52l [102L if 
£ m M = M £ m . We proceed to extend this concept to semi-infinite matrices 

Definition 2.12. A block semi-infinite matrix M is persymmetric if 

r|M = M T r|; 


i.e., if its blocks satisfy £ m M mi] = (M [l])[k] ) T %]|. 

Notice that in the diagonal, k = l, one recovers the standard persymmetry property of the diagonal 
square blocks. 

Proposition 2.12. The moment matrix G is a persymmetric semi-infinite matrix 

(2.9) rjGr] = G T . 


Proof. To prove it we perform the following sequence of equalities 


rjGri=o r)x(z(0))dM.(0)(x(z(0))) f ri 
Jt d 

= o x(z(-0))dq.(0)(x(z(-0))) 

Jt d 


t 


= o x( z (9)) d p(0)(x(z(0))) 
Jt d 


T 


= (o x( z (d)) d p(0)(x(z(0))) 
v Jt d 


T 


=g t . 


□ 


Proposition 2.13. The following properties hold true 

(1) The matrices H are persymmetric r|Hr| = H T and 

(2.10) § = r|Sr|. 

For real measures we have r|Hr| = H. 

(2) The following is satisfied 

r||3r| = 0. 

(3) The MVOLPUTfulfill 

(2.11) r|T>(z) = ), T)^(z) = T>(z _1 ). 

Proof. From the persymmetry property ( |2.9| we get 

qG =G t ti, 

and using the Gauss-Borel factorization ( |2.5| we obtain 

gpS -1 = H T (S T ) - 1 Ti$ t l-r 1 =q, 

and, as we have a lower triangular matrix on the LHS and a upper triangular matrix on the RHS, the only 
option for is to be a diagonal matrix, i.e., equal to the reversal matrix r\. □ 

We discuss now some interesting matrices, that we name as partial exchange or reversal matrices which 
are useful in the finding of interesting parity properties of the MVOLPUT. For that aim we need to intro¬ 
duce 

Definition 2.13. We consider the signature matrices I Q e C DxD , a £ (1,..., D), these are diagonal matrices zvith 
their diagonal coefficients being 1 but for the a-th entry which is —1. 

With the help of these signature matrices we define 
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Definition 2.14. We consider partial reversal matrices r) a , a G {1,..., D} which are block diagonal semi-infinite 
matrices with coefficients given by 




^a| k) ,I a ai k) ^I a a| k ',a[ k ’ 1 ' 


In Appendix |B] we give some examples for the cases D = 1,2 see . 
Proposition 2.14. For a, b e {1,..., D } the partial reversal matrices fulfill 


hahb = hbha, 

Proof It is easy to realize that 

I Ml 


hi = i, 


=r i- 




r=1 


(k) (k) 

*• ,a r 


(hbl 


,(k) „(k) 

A-r 


= 5 


Iuaj kl ,Iba[ k * ''I a IbO£j' iJ ,a{ 


= 5. 


Q=1 


„(k) Hk). 


The signature matrices I a are square roots of the identity 1^=1 and so are the corresponding partial 
exchange matrices r| a . Finally, from 


I Ml 

L 


„(k) 




n 


( T l 2 )„(k) (k) • • • (ho) „(k) fk) - 5 IK.J 

«t, t a r 2 a r D -)> a l *1 lD«j i a 


r l v) r D-1 =1 

but h • • • I D = —I and therefore we recover the complete exchange matrix q. 


(k) „(k), 
■j 


□ 


Definition 2.15. (1) We identify {1,..., D} = Zq as the D cyclic group Zd, being D the identity for the sum 

a + D = a (mod D). 

(2) For each variable a G Z ^ we introduce the partial reversal operator V a : (C*) D —> (C*) D given by 

a(^l i • • • i Zat • • • > ) - = (^-l) • • • > ^.q ) • • • > ZD ) ■ 

(3) This partial reversal operator extends to the algebraic torus (C* ) D the partial parity operator in the unit torus 

P Q : T° -a T d , (0i,...,0 a , ...,0 d ) •—> (0i,...,—0 Q , ...,0 d)- 

(4) Its action on the measure will be denoted by V* d p(0) = d \i{V a Q ). 

(5) Given a subset cr g 2 Zd we consider the reversal operator V a := Pa defined in the algebraic torus, the 

aG a 

partial reversal matrix q a := q Q , the complex vector z ff := P ff z G (C*) D and the corresponding moment 

aCc 

matrix := G -p*^. 

(6) The LU. factorization G (<T -' = ^ leads to the corresponding partial reversal trans¬ 
formed MVOLPUT 0 (cr) := S^x an d := §^x an d second kind functions := (^(S (o '*) ] ^x an d 

£{*) := 

(7) For a subset a g 2 Zd we consider its complement Ccr defined by the following two conditions: aljCa = 
(1,2,..., D} and cr p| Co = 0. 


Proposition 2.15. (1) Given a subset cr C Zd the partial reversal matrices q ff act on the monomial vector as 

follows 

haX(z) =X(Zff)- 

The above property models, in unit torus variables 0, the partial parity operation of reversing the signs of 
appropriate angles 0's. 

(2) The follozving equations holds true 

G [ff) =q a Gq ff =q Ca G T q Cff = (G ,CtT) ) . 


16 


GERARDO ARIZNABARRETA AND MANUEL MANAS 


(3) The partial reversal transformed MVOLPUT and second kind functions are connected to the original MVOLPUT 
and second kind functions by 

=r| C(T &(z C(J ), =q CtT <f)(z c J, C [a] {z)= qc^ZCff), ^{z) =^] C(y C[z Ca ). 

Proof. We only prove the third statement being the others straightforward. Assuming the factorization 
G* 0 '* = (S^) M and recalling rj^GriQg. = (G ,ff *) T we conclude 


= § (<j) nn ‘T 1 = 


hCcr^ hCo"- 

Therefore, 

s (<J) =ncff§nc<n 

and the result follows. 


(h ( ')) t (flicks 1 " 1 )"’ 


T 


H . 


h 11 =h C a H nca) 


^ (<I) =hCaShCc 


□ 


2.6. Plemej type formulae for MVOLPUT. Here we discuss integral expressions for the second kind func¬ 
tions and Plemej type formulae. We follow the papers by Mohammed H72l and If73] . 

Definition 2.16. For each subset a G 2 Z ° we define its right boundary by 

3a :={i G Zd : i G a and i + 1 G Ca}. 

D 

Then, we introduce the orthant polydisk as the polydomain (D d )o- := X At with A; := C\D/fi G a and A t := D 

i=1 

when i G Ca and define, for a f Zd, the corresponding integer orthants 


(Z D ) ff :=XZ|, 

i=l 


D 


Zi := 


Z-, 

i G (a\ 3a), 

z<, 

i G 3a, 

Z+, 

i G (Ca\ 3(Ca)), 

z>, 

i G 3Ca, 


zvhere Z^ := Z± \ {0}. For a = Zq we define 

(Z d ) Zd := (Z_)° \ {0}, 

where 0 G Z D is the zero vector. 

Proposition 2.16. The set of multi-indices splits into disjoint integer orthants 

Z D = [J (Z D ) ff , with (Z D )o- n (Z D ) (J / = 0/or a f a'. 

ce2 z D 

Proof. See Appendix |C.3| 

Definition 2.17. For each subset a G 2 Z ° we define the Fourier orthant vector 

_ J 0, «GZ d \(Z%, 

\z“, 

and the Cauchy-Mohammed kernel C a {z, C),for £ g T D and z G (D d )o-, asfollozvs 


□ 


Proposition 2.17. We have 


Cff(z, C) :=(x(0) t Xa(z)- 
X= Y- Xa ' 


ce2 z D 


Proof. It is a consequence of Proposition|2.16| 


□ 
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Proposition 2.18. The Cauchy-Mohammed kernel has the following expression 


C<r(z, C) = ^ 


n 


Zi 


_ J ie(cr\3(j) 


Zi Ci 


Zi 


n z,-C- 

i£(o\3o') Za Li 


n 

ig3cr Zi Li 

-1 


n —- 

ie(Ccr\aCcr) & — Z i 


n 


Zi 


igSCcr & Z i- 


CT / Z D , 
O = Zd, 


/or C e T d and z E (D D ) ff . 


Proof. For a / Zq from the Definition 2.17 we have 

Ccr(z, C) = X r “ z “ 

ae(Z D }a 

=nt l gxl 


i=l aiSZi 


Now, we deduce that 


GXJ 

£ (Ciz/ 1 ) 1 / iE(ff\3ff), 


n=0 

oo 


L (CiZi 1 ) n , iG3a, 

L f-CX-i ~ a i J Tl=1 

'H N oo 

£ (C/ 1 ^) 11 , i E (Co \ 9 (Cct)), 


aiSZi 


n=0 

OO 


L (C t Zi) n , iESCct, 


n= 1 


hence, using that z e (® u )a and the basic results on geometrical series we find 


C ff (z, C) = ( Yl 


1 


1 - Ciz: 


-1 


n 


CjZ t 1 
1 - CiZ/ 1 


n 


i 




n 


ti’zi 




□ 


iS(cr\9(j) ig3(j '" a X ig(C(r\3Coj ^ 1 igBCff 

and we get the result. For cr = Zq the result follows similarly. 

This result shows that these integral kernels (—1 )^C a coincide with the ones in 172. 73 1. 

Definition 2.18. Given z E (D d ) ct we introduce the corresponding orthant components of the second kind functions 

C a [z) :=(S -1 ) f x o-(z), &(z) ^(S'Yxafz). 

Observe that we could also write 


Proposition 2.19. For z E 




the orthant components of the second kind functions are integral transforms 


according to the Cauchy-Mohammed kernels 

C a [z) =(H _1 ) f o ^(e 10 )dfi(0)Cff(z,e ie ), 

Jt d 


^ct(z) =H 1 o O(e 10 )dp(0)C (7 (z,e 10 ). 

Jt d 


Proof. It follows from Definition 2.18 that 


Co-(z) =(H 1 ) t cj) §x(e 10 )dp.(0)(x(e 10 )) t Xa(z), &(z) =H 1 <b Sx(e 1B ) d p(0)(x(e 10 ) ) Xc[z). 

J T D V / J T D 


-i 


je 


t 


□ 


and Proposition |2. 1 8| gives the desired result. 

Definition 2.19. Given a function f : T° —> C the function denotes its holomorphic extension to (D D ) a having 
f as its boundary value at T°. 
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Theorem 2.1 (Plemej formulae). Take an absolutely continuous measure dp(0) = (1(0) d0 and assume that 
Holder continuity for O a jl, 6 a (l e C a (T D , C), 0 < a < 1, holds. Then, we have 


ze(D D ),T 

z->£ 


lim cb $ a (e I0 )C ff (z,e lb )d(l(0)+ lim cb 


J 0 


t d 


ze(D D ) Cff 
z— 


J0\ 


j e 


t d 


= (2 7 t) D ((a a A) w (o + (^aA) [Cff] (c)), 


and 


z£ 


lim cb O a (e 10 )Ca-(z,e 1B )dM.(0) + lim cb iDJe'^CcJz^^ld^©) 


J e 


T D 


z->£ 


ze(D D ) Cff 
z— 


_i 0 'i 


10i 


T D 


= (2 7 r) D ((d> ot p.) [CT] (C) + (O a bf) [Ca] (0), 


for C e T u and a € Z . Equivalently, 


Y lim o ^ a (e 10 )C lr (2,e lfl )di‘(0] =(27t) D d a (C)|l(C), 


Y Hm o ®«(e le )C„(z,e ie )d H (0) =(27t) D ® a Mfl(C). 


Proof. For the proof we need some observations regarding ||72| and 11731 . Obviously, the MVOLPUT belong 
to the Wiener algebra of the D-dimensional unit torus, 

®a, &a € W D = |f(z) = Y_ c « z “ : ||f||w D : = Y- l C «l < °°}- 

aez D aez D 


We assume that dp = (1(0) d 0 with |1 e W D , and consequently, given that we are dealing with Laurent 
polynomials, ® a |2, ® a p £ W D . Following l72l[73l for any Holder continuous complex function on the unit 
torus cp e C a (T D ,C), 0 < a < 1, one defines its integral transformations ([l' 11 and its boundary integral 
conjugates cbl 0 ' as follows 


4) (tT) (z) := 


(27t)t 


■ o 


T D 


cp(e 10 )C ff (z,e 10 ) d0, 


1 


4>* (n) : =-d ° ^(e 10 ) 0 ^ 10 ’ 1 !) d0, 

7T J-jD 


where 


C(t,n) = bl t 


“Ci-m 


Observe that 


4> ftT) C“n) =io 4jl ff, (e 10 )C(e 10 ,T|) d0, 

71 J'j’D 


and for £ e T u we have Plemej's type formulae 


4) (ff) (0 + 4> (CtT) (0 =cp [cT] (C) + tp lC(l] (C), 
X. ^i a) (C)=2 D 4)(C). 


ae2 z D 


ze (P D )a, 
q GT d , 


71 E T L 


where the holomorphic extension <p ‘ r (z) to (B D ) (r is obtained from the Fourier series of cp by disregarding 
the mult-indices not in (Z D ) ff ; this holomorphic function in (B D ) ff has cp as its boundary value at the unit 
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torus T d . Thus, the orthant expansion 

ch (CT) (C) = 2 D cp(C) 

cf£2 z d 

holds. □ 

This Plemej formulae can be recasted for the second kind functions as follows 
Corollary 2.1. For any point in the unit torus, £ £ T°, the MVOLPUT and its second kind functions satisfy 
Y_ Hm C a (z) =(27t) D Hta(£)ftS), lim £ a [z) =(27r) D HO(C)WC), 


ze 

a£2 z D z ^ c 




or equivalently 


lim C ff (z)+ lim C c<T (z)=(2 7 r) D Ht((aA) W (C) + (aA) [Cff] (C)), 
e(B D ) CT ze(B D ) Cff v J 

z—z— 

lim £ a [z)+ lim d c<T (z)=(2 7 t) D H((d)A) [ll] (C) + (OA) [Cl7] (C)). 

P(I) D L ZP(l) D )n_ V / 


z—>£ 


zeOD> LJ ) CCT 
z— 


2.7. String equations for the moment matrix. 

Definition 2.20. We introduce the spectral matrices T a , a £ {1,..., D}, which are tri-diagonal block matrices 

0 ...\ 


with entries given by 



f 0 

(Ta) [0],[1] 

0 


( T a)[1],[0] 

0 

( T a)[1],[2] 

— 

0 

( T a) [2],[1] 

0 


0 

0 

(Ta) [3],[1] 


V i 



(Y a ) to „(k+n 

i 

:=6 (kj, 

a) +e a 

,aj k+1) = 

(T a )„(k+n (ki 

: - 8 «' k+1 ) + e Q ,«; k) - 


'■J 




> a (^i) )a w_ eQ . 


where {e Q }° =1 is the canonical basis in C D . 

Proposition 2.20. The spectral matrices satisfy 

T 'Y' r \r 'Y' 'V'—1 v'T „ 'Y' 

a ' b — I b * a» 'a — * a > % * b — 

and the very important spectral property 

T a x(z) =ZaX(z )• 

Proof. See Appendix |C.4[ 

For examples with D = 1,2 see Appendix [B] 

Definition 2.21. For any a = (<xi,..., ocq) £ TP we introduce the matrix 


T b ri a a/b, a,b £ {1, — , D}, 
YkV a = b, 


□ 


( 2 . 12 ) 


Ta :=fj(T a ) aa . 


a=1 


Proposition 2.21. For the spectral matrices the following properties hold 

(1) Given a multi-index a £ Z D we have T_ a = (T a ) _1 = (Y a ) T = r|T a ri; i.e., Y a is an orthogonal 
persymmetric block matrix. 
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(2) Y a is a 2\a\ + 1 -diagonal block matrix, with at most |a| block superdiagonals and |a| block subdiagonal. If 
|a| is odd, the main diagonal and the even superdiagonals and subdiagonals vanish, if\a\ is even then the odd 
block superdiagonals and subdiagonals are zero. 

(3) For any couple of multi-indices cq, ct 2 £ we have Y a , +a2 = T a , T a2 . 

(4) The following eigen-value property holds T ax{z) = z“x(z). 

Proof. (1) From ( |2.12[ ) and the properties of the exchange matrix r\ we get 

D 

= n (h^ah)“ U 

a=1 

and recalling that {Y a }^ =1 is an Abelian set of real persymmetric orthogonal matrices; i.e., r|Y Q r| = 
Y (1 = Y~ 1 we deduce that we get 

nr a 11 =n(r:r=n m - * 0 =n ( T ar aa 

a=1 a=1 a=1 

_Y'T at—1 y' 

~ ' a ~ 1 a — 1 - oc • 

(2) As Y a are tri-diagonal block matrices then the product of |a| of them will be a 2| a| +1 -diagonal block 
matrix, and recalling that the main diagonal vanishes in Y Q we deduce the rest of the statement. 

(3) It is a consequence of ( |2.12| and the commutativity of the Y a among them. 

(4) Follows also from the eigen-value property described in Proposition |Z20| 

□ 


We are now ready to show a very important symmetry 
Proposition 2.22 (The string equations). For any multi-index «eZ D the moment matrix satisfies 
(2.13) Y a G = GY a . 


Proof. It is a consequence of ( |3.3[ ) and the following sequence of equalities 

Y a G =Y a cj) x( z (9)) d p(0) (x(z(—0))) T 


t d 

= o Y a x(z(0)) d p(0)(x(z(—0))) T 

Jt d 

= o e ia ' e x(z(0)) d p(0) (x(z(—0))) 1 

Jt d 

= o x( z (S)) dq(0)(e ia ' e x(z(-0))) 

Jt d 

= o x(z(0)) dq(0](Y^x(z(-0))) T 

Jt d 

= o x(z(0)) dq(0)(x(z(-0))) T Y a 
Jt d 

=GY a . 


2.8. Three-term relations. The symmetries of the moment matrix suggest 

Definition 2.22. We consider the Jacobi matrices 

J«:=SY a S"\ T«:=§Y«$^ 

and the reversal Jacobi matrices 

C a :=%T a S-\ 

Proposition 2.23. The following properties hold true 


□ 
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( 1 ) = 

(2) J a is a [2\oc\ + 1 )-diagonal block semi-infinite matrix 

(3) We have Jo = I • 

(4) For any couple of multi-indices ai,a 2 G Z D we have J a ,+ a2 = Ja, Ja 2 - particular, (Ja)^ 1 = J-a 
V«GZ°. 

(5) We /ifloe 

C« = II C a H Co =tl, = Sr|Y a (§) 1 , 

and consequently C a is a (2|a| + 1 )-diagonal block semi-infinite matrix. 

(6) The following properties are satisfied 

llja = c a , L] Ja = C a , LlJaT) = J— a . 

Proof. See Appendix |C.5[ □ 

Proposition 2.24. The Jacobi matrices satisfy the following eigen-value relations 

J a 0(z) = z a O(z), T«6(Z) = z«d(z), 

and the following relations are fulfilled as well 

C a O(z) =z a $(z- 1 ), C a 1 $(z) = z a ®(z -1 ). 


The basic Jacobi matrices J a = SY a S 1 are tri-diagonal block semi-infinite matrices 


Ja 


/ (Ja) [0],[0] (Ja)[0],[1] 0 

(Ja)[1],[0] (Ja)[1],[1] (J a) [1 ], [2] 

0 (Ja)[2],[1] (Ja)[2],[2] 

0 0 (J a) [3], [2] 

V : : 


0 

0 

(Ja)[2],[3] 
(J a) [3], [3] 


0 ...\ 

0 

0 

(Ja) [3],[4] 


with blocks given in terms of quasi-tau matrices and first subdiagonal coefficients by 


(2.14) 


(J a) [0] [0] --(^a)[0],[1]P[1] =- H [O]0[l ] C^ a )[1],[ O ]H [o] 1 , 

(Ja)[k],[k—1] = H [k] (a) [k], [k— 1 ] H [k —1 ]) 

S (J a) [k] [k] = p[k](^a)[k-1],[k] ~ (^a)[k],[k+1] P[k+1] 

= H[k] [(YaJ^^k-ljP^ — |3q c+1 ] (^a)[lc-Ll],[lc]] Hj-^ 
, (Ja) [k],[k+1] = C^a)[k],[k+1] 


Similarly, J Q has the same form and the blocks are obtained from the previous formulae with the replace¬ 
ments 


H <—¥ H f , 


(3 <—» (3. 


The matrices C Q are tri-diagonal with blocks given by 

(Ca)[0][0] =-(b^a)[0]mP[1] = —H-|p] P|^] (TITa)[1 ][0]^[0] 1 > 

(Ca)[k],[k-1] = l~l[k] ("n^ci)[k],[k—— 1 ]? 

(C Q )[k][k] = ^[k]( 1 lla)[k-1],[k] _ ( r l^'a)[k],[k+1]P[k+l] 

= H [k] [(hTa^kUk-llpM - Pj+u^YaJlk+llJkllllpcp 

(C a )[k],[k+1] = h^a)[k],[k+1]) 

and Cq 1 is also tri-diagonal with blocks as above with the replacement given by the following prescription 


H <—¥ H t , 


(3 <—> (3. 
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Definition 2.23. With the notation 

T-a-.= T-\ a G {1,..., D}, 

the vector n = (n^i,..., n D , no,..., ni ) T 6 C 2D and ft := SznP we consider: 

(1) associated longitude one homogeneous Laurent polynomials 


D 


L n (z) := Y_ —“ + Y naZa ’ 

i Z a i 

a=1 a=l 

(2) associated slashed semi-infinite matrices 

D 

ft = Ln(T) := ^(n-aT- 1 + n Q T Q ), 


a=l 


D 


a=1 


Lfl(z) := Y_ n -a z a + Y — • 


D 


ft = Lft(T) := Y + Y < T a' • 

a=1 a=1 a=l 

Observe that 

Proposition 2.25. The following formulae holds 

ft ] =ft, ft x(z) =n(z)x(z). 

Proposition 2.26 (Three-term relations). The MVOLPUT satisfy the following recursion type relations 
hnM^Ik] ( z ) = H[ k ] ft [k] ) [k_i]H^_i ] ())[ k _ 1 ](z) + (|3[k] ft [1c—1],[lc]— ft [k],[k-1]P[k+1])4 ) [k]( z ) + ft [k],[k+1]4>[k+1] ( z )> 
hn(z)$ [k ](z) = (H [k] ) t fi [k],[k— 1 ](hT [lc 1 _ 1] ) t $[i c _ 1] (z) + (j§[k] ft [k-l],[k]- ft [k],[k-1]fWl])$[k]( z )+ ft [k],[k+ 1 ]$[k+ 1 ]( z )> 
and its conjugate partners 

hn(z)$ [k] (z _1 ) = HjjOl ?f)[k] > [k- 1 ]H[ k L 1] ct>[k- 1 ]( z ) + (i[k]h fl)[k- 1 ],[k] - (T1 l%],[k+1]P[k+1])4>[k](z) 

+ (h 1^)[k],[k+1]4 ) [k+l]( z )> 

hn(z)4>[k] ( z_ 1 ) = H[k] (ri 7i)[k],[k-1](H[k-l]) 1 $[k—1] ( z ) + (P[k](Tl ft) [k—1],[k] ~ (h ft) [k],[k+1] P[k+1]) $[k] ( z ) 

+ Ol T^-) [k],[k+1]‘t’tk+l] ( z )- 

2.9. Christoffel-Darboux formulae. 


Definition 2.24. The Christoffel-Darboux kernel is defined by 


k—1 

K (k) (zi,z 2 ) := Y ($[i]( z i)) f hT^] 1 [i]( z 2 ) 


1=0 


Proposition 2.27. Both sets of MVOLPUT {cf)q } | ct |< k and {$q}| a |< k are linear basis for the linear space QJz^ 1 ] 
of Laurent polynomials of longitude less or equal to k. 

Thus, a Laurent polynomial L can be always expanded as 

L(z) = Y 

aeZ D 

0^|a|<Coo 

= y~ '^a4 ) a( z )> 

aez D 

0^|a|<Coo 


with coefficients % a , A a 6 C. 

The corresponding projections are 

Definition 2.25. We introduce the two projectors onto [z ±] ] as the following k -th truncated sums 
rr (k) fL):= A«4>«(z), fl< k )(L):= Y X«$«(z), 

aez° a€Z° 

0^|aj<k 0^|a|^k 
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Due to the biorthogonality \2.7) we have 

Proposition 2.28. (1) The coefficients in the MVOLPUT expansions of a Laurent polynomial L are given by 

L(e 10 )4> a '(e ie ) d p.(0) 


^oc — 


Y_ L(e ie )$ a ,(e i0 )(H- 1 ) a , a dp(0), (FT 1 )*,*' 


|a'|=|a| 


Jt d 


|a'[=|a| 


T d 


(2) The projections given in Definition 2.25 can be expressed in terms of the Christoffel-Darboux kernel as 


TT (k) (L)(z)=o L(e i0 )K (k+1) (e i0 ,z)dp(0), fl (k) (L)(z) =o K( k + 1 )( z , e i0)L(e i0 ) d (t(0). 

Jt d Jt d 

Therefore, these kernels are the integral kernels of the projection operators in the space of multivariate 
Laurent polynomials of longitude equal or less than k. This interpretation leads at once to 

Proposition 2.29 (Reproducing property). K^jzi^) = J t d K (k) (zi,e 10 )K( k ^(e 10 ,Z 2 ) d |x(0). 

From the Gauss-Borel factorization we deduce 

Proposition 2.30 (ABC Theorem). An Aitken-Berg-Coller type formula 

K«(z 1 ,z 2 ) = ( X [k] (zi)) t (G [k] r 1 X [k] (z 2 ), 

holds. 

Proposition 2.31. The Christoffel-Darboux kernel satisfies 

K (k) (z 1 ,z 2 ) = K 

Proof. It is a consequence of ( 2.\1\ and the per symmetry of H: 

k-1 

K ,k) (z!,z 2 ) = Y_ (43[i](zf 1 )) T «F| [l ]|H [l] 1 (j) [l] (z2) 


1=0 

k—1 


($ [l] (Z 1 ))' t H jy 1 £j [y I $ [y (3 


.-1 


=0 




□ 


Proposition 2.32. The Christoffel-Darboux type formulae hold 
K [k] (Z!,Z2) = 


[k] ^ ($[k](zi)) t yi [k],[k-i]H [k 1 _ 1] c|)[k_i](z 2 ) — ($[k-l](zi)) t H [k 1 _ 1] yt [k-i],[k]4 ) [k](z2) 


L n (z 2 ) - L n [z^) 

l^)[k],[k—l]hl^J ^[k—l](z 2 ) - (4?[k^l] Cz^ 1 )) T (H^ 1] )- 1 (q 7f)[ic i] >M 4>[k](Zl) 


i-1 ’ 


f-n(Z2) L n (Z-| ) 

_ ($[k](zi)) t (Tfq)[k],[k- 1 ](Hj_ 1] ) _1 $ [ k-i ] (zJ 1 ) - ($[k— 1 ] (Zi )) t H [ j c ^y [yi q)[k—l],[k]$[k] 1 ) 

L n (z 2 ) - Lnjz^ 1 ) 

Proof. In order to get the first of the three expressions one has to consider the two possible ways of com¬ 
puting the expression (cf) [k+1] (zi))t [(11 -1 (L n (J)) k ' 1 ] <J) lc 111 (Z 2 ) letting the operator between brackets act 
on the polynomial located on its right or left hand side; i.e., 

($ [k+1] (z!))t [(H- 1 (L n (J)) [k+1] 4> [k+1] (z 2 )] =L n (z 2 )($ [k+1] (z 1 ))t( H [k+1] )- 1 4> [k+1] (z2) 

~ ($ [k] (■Z1)) t H ('T a) [ k ] [k+ ! ] (j) [k+ ! ] (Z2), 

'($ [k+1] (z 1 ))t(H- 1 L n (J)) [k+1 l] 4) [k+1] (z 2 ) = L n (z 1 - 1 )($ [k+1] (z 1 ))t(H [k+1] )- 1 4) [k + 1] (z2) 

— ($[k+1] (zi)) t (T a ) [ k+1 ] [lc ] H^l 4)[k] (Z 2 ) . 
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Since both results must coincide, subtracting both we get the desired expression. The remaining two equal¬ 
ities can be proven in two different ways. The first and simpler one is to use the expression we have just 
proven and remember the relations between <[) and $ established by q. A different approach is to compare 
the two possible ways of computing 

(4> [k+1] (zf 1 C Q ) [k+1] }4) [k+1] (z 2 ) 

and 

($ [k+1] (zi))t 

□ 


([C- 1 ] T [H- 1 ] T ) 0t+1] 4> L * +IJ (Z2). 


. k I 1 / 


3. Darboux transformations 

We now discuss perturbations of the measure given by the multiplication by a Laurent polynomial. This 
is the extension, in the multivariate Laurent orthogonal polynomial setting of the Christoffel transforma¬ 
tion, or in a geometric-differenttial setting of a Levy transformation, or in the theory of integrable systems 
of a direct Darboux transformation. 


3.1. Laurent polynomial perturbations of the measure. 

Definition 3.1. Given a multivariate Laurent polynomial the corresponding Darboux transformation of the measure 
is the following perturbed measure 

(3.1) Td |x(0) = L(z(0)) d |x(0). 


Proposition 3.1. The perturbed measure L d p is real if and only if L[z 1 ) = L(z); i.e., when the Laurent polynomial 
can be written as 


l I M|/2 

L(z) = U + £_ X (L 

k=1 i=1 




,(k)Z 


+ L (ic)Z 

a. 




When we have 


l I Ml/2 


l»> 2 L L |L „™i 

k=l i=1 


definite positiveness of T d p is ensured. 


Proof. The perturbed the measure can be written 

l I Ml/2 

L(z(0))dp(0) = (ip + 2 Y Y_ |L a (k)|cos (a{ k) ■ 0 + argL a(k ))) dp(0). 

k=1 i=1 1 1 

from where the sufficient condition for positiveness follows. □ 

Another manner of ensuring positivity is as follows. For D = 1 the Fejer-Riesz factorization [39/ 77] 
allows us for expressing a Laurent polynomial which is non negative in the circle T in the following manner 
L(z) = Q(z -1 )Q(z), for a polynomial Q; notice that this expression when evaluated on the circle T takes 
the form L(e 10 ) = |Q(e 10 )| 2 . In the multivariate scenario the situation is quite different. As proven in 
Il36l| and discussed further in Il48l Fejer-Riesz factorization can be extended; indeed, we can write any 
multivariate Laurent polynomial strictly positive in T° as the Dirtschel's finite sum of squared magnitudes 
of multivariate polynomials 

r r 

(3.2) L(z)=^Q i (z- 1 )Q i (z), L(e ie ) =^|Q i (e 10 )| 2 . 

i=1 i=1 

Proposition 3.2. The moment matrices TG and G satisfy 

TG = L(T)G = G L(T). 
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Proof. We prove the first equality 
(3-3) 


T 


(3-4) 

(3-5) 

(3-6) 


TG =o x(z(e))(Tdp(e))x(z(—0)) 

Jt d 

= 1 x( z (9))(h(z(0)) d |x(0))x(z(—0))" 

Jt d 


= o (L(T) X )(z(0)) d |t(0)x(z(— 0)) 
Jt d 

=L(T)G. 


T 


from Definition 12.61 
from ( |3.1| 
from ( 2.20} 


The second equation follows similarly. 

We proceed to introduce a semi-infinite matrix that models the Darboux transformation 
Definition 3.2. The resolvents are 

tu :=(TS)L(T)S _1 , d>t := (§- 1 ) t L(T)(T§)t. 

We also introduce 

Definition 3.3. The adjoint resolvents are 

A:=S M:=STS~\ 

Proposition 3.3. The resolvents and the adjoint resolvents satisfy 

=(TH -1 )cuH, M =Hd5 t (TH) _1 . 

Proof Proposition |3.3| and the corresponding Gauss-Borel factorizations lead to 

(TS)- 1 (TH)((T§)- 1 ) t = L(r)S -1 H(§ -1 )t = S^H^^KT), 
and an appropriate cleaning of these equations to 

(TH)(§(T§r 1 ) t = (TS)L(Y)S _1 H, S(TS) _1 (TH) = H(§ _1 ) t L(T)(T§) t , 

from where the result follows. 


□ 


□ 


Proposition 3.4. The resolvents cu and d) are block upper triangular matrices having all their block superdiagonals 
above the m-th block superdiagonal equal to zero. The adjoint resolvents M and K/l are block lower unitriangular 
matrices having all their block subdiagonals below the m-th block subdiagonal equal to zero. 


Proof. The reader should notice that, if the longitude of the perturbing Laurent polynomial is l (L) = m, the 
semi-infinite matrix L(T) is a (2m T- 1)-banded matrix with only the first m block superdiagonals, the main 
diagonal and the first m subdiagonals different from zero; hence. Definition |3.2| implies that cu has all its 
superdiagonals above the m-th superdiagonal equal to zero and that dX has all its subdiagonals below the 
m-th subdiagonal equal to zero. Observe also that Definition |T3| implies that the adjoint resolvents M and 
M are block lower unitriangular matrices. Now, Proposition 
gives the satated result. 


3.3 


linking resolvents and adjoint resolvents 

□ 


Regarding the deformed MVOLPUT and the original ones we have 


Proposition 3.5. The following relations among Darboux perturbed and original MVOLPUT hold true 

(3.7) cuO(z) = L(z)T<D(z), d)d(z) = £(z -1 )Td(z), 

(3.8) MTO(z) = 0>(z), iQlTd(z) = 6(z), 

Proof. To prove equations ( |3.7| ) just observe that 

cu4>(z) = (TS)L(T)S _1 S X = (TS)L(T) X (z) = L(z)T®(z), 
d>a(z) = (TSXLmjt^Sx = (T§)(L(T))t x (z) = L(z _1 )t6(z). 

The relations ( |3.8| > are a consequence of Proposition |33] □ 









26 


GERARDO ARIZNABARRETA AND MANUEL MANAS 


3.2. Interpolation, sample matrices and poised sets. 

Definition 3.4. The algebraic hypersurface of zeroes of the Laurent polynomial L in the algebraic torus is denoted by 

Z(L) :={p<E(C*) D :L(p)=0}. 

Definition 3.5. A set of nodes in the algebraic torus 

A4,m :={Pj}£f C(C*) D 

is a set with rk )1Tl := Nk+ m -i — Nk_i = |[k]| + • • • + |[k + m — l]|nodes. Given these nodes we consider the 
corresponding sample matrices 


(3.9) 


I 


m ._ 

k •— 


V^fc+m-llCPl) 


^[k] (Pr^mJ 
‘t’fk+m—1] (Pr k>m ) 


G ■ x ' 1 '' c ’ rrL 


I[k,m] :=(4)[k+m](Pl),---,43[k+m](Pr k , m )) G C |[1c+m]|xrk --. 


Lemma 3.1. When the set of nodes A/k,m C Z(L) belongs to the algebraic hypersurface of the Laurent polynomial L 
of longitude m = £(L), the resolvent coefficients satisfy 


^[k],[k+m]^-[k,ra] ~h ([k],[k]> • • • > ^[k],[k+m— 1])^-k 9. 

Now we introduce an important type of node sets within the algebraic hypersurface of the Laurent 
polynomial L of longitude m = £(L) 

Definition 3.6. We say that the set A/k,m C Z(L) is a poised set if the sample matrix given in ( |3.9| is non singular 

det IJJ 1 f 0. 

Theorem 3.1 (A Christoffel formula). For a poised set of nodes A/"k,m C Z(L) the Darboux transformation of the 
MVOLPUT can be expressed in terms of the original ones as the following last quasi-determinantal expression 



f 

4>[k](z) 

ti (L(Y)[k] [k+m] 0 

T4)[k](z)= ©* 


4>[k+m~1]( z ) 


V Z[k,m] 

^[k+m] (^0 


Proof. Observe that Lemma 3.1 together with eu[ k ]ji c+m ] = (LfYj) , , implies 


([k],[k] > • • • ) ^[k],[k+m—1]) — (L(Y)) [^[k+ml^dk™-] (^-k ) 


and from L(z)Tc|)(z) = cuc[)(z) the result follows. 

Proposition 3.6. The resolvent coefficients can be expressed as last quasi-determinants 



/ 

°[k],[k+l] 



0[k+l-1],[k+l] 

^[k],[k+l] = (I-('V"))[ k> ] j [| c _ HTn j©* 


I [k+l] 



0[k+l+1],[k+l] 


V Z[k, m ] 

®[k+m],[k+l] 


for l = 0,..., m — 1. 


□ 
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Proposition 3.7. We have the following Darboux transformation formulae for the quasi-tau matrices 


( 


Hfc] \ 

°[k+1],[k] 


TH [k] =(L(T)) [k] 

,[k+m]®i 

and for the transformed first subdiagonal coefficients 

(3.10) (T(3) [k] (L(T)) [k _ 1] 

,[k+m— 1] = (L(T)) [k],[k+m— 1] + L(T)) 

[k],[k+m]®* 

Observe that ( |3.10| can be written 

(TP)m = (h(T)) [k];[k+m _ 1] (L(T))+ +m ^ 1])[k _ 1] 


V I [k,m] 

°[k+m],[k] J 




/ 

0[k],[k+m—1] 

+ L(Y))[ k ] )[k+m ]@* 


0 [k+m—2], [k+m+1 ] 
^[k+m+1] 



V I [k,m] 

P [k+m] 

( 

°[k],[k+m-1] 

\ 


+ L(Y)) [k ] )[k+m] 0* 


ym 

A k 


V ^-[k,n 


0 [k+m—2], [k+m+1 ] 
^[k+m+1] 


P [k+m] 


(L(T)) + 


[k+m— 1 ], [k— 1]* 


whenever the right inverse (L(Y))jj c+Ta ^ - k 1 of (h(Y)) [k _ 1])[k+m _i] 


exists. 


Proposition 3.8. The evaluation of the Laurent polynomial L on the Jacobi matrices J = (Ji,..., Jd) t or on its 
perturbations TJ = (TJi,... , TJd) t are linked to the typical Darboux LU and UL factorizations of the Jacobi and 
Darboux transformed Jacobi matrices ; i.e., 


L(J) = Mcu, 


L(TJ) = ojM, 


respectively. 


Proof. From Definition 3.2 we get cu = (TS)S _1 SL(Y)S _1 = (TS)L(Y)(TS) _1 (TS)S _1 and therefore, using 
Proposition |3.3| and Definition 2.22 we get Mcu = L(J), from the first equality, and tuM = L(TJ) from the 
second equality. □ 

From the first equation in the previous Proposition we get 

Proposition 3.9. The block truncations (Q(J))^ admit a LU factorization 

(L (J)) [1c] = M [k V k] 

in terms of the corresponding truncations of the adjoint resolvent and resolvent tuH 
Proposition 3.10. The truncated matrix (L(J))^ is regular and 


k—1 


det ((L(J)) [k] ) = ]7] 


detTH 


[i] 


1=0 


detH 


[i] 


Proof. To prove this result just use Propositions |3.8[ the explicit expression 

= ( TH [k])H [k] 1 , 

and the assumption that the minors of the moment matrix and the perturbed moment matrix are not zero. 

□ 
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3.3. Nice Laurent polynomials. Flavors from tropical geometry. So far, for the analysis of Darboux trans¬ 
formations of MVOLPUT, we have followed the approach given in Ill4l for MVOPR. However, for the 
polynomial ring C[zi,... ,zq] we had an ordering, the graded lexicographical order, wich is also a grad¬ 
ing: deg(PQ) = degP + deg Q which is missing in this contex t. N ow, for the Laurent polynomial ring 
C[z^ 1 ,..., Zq 1 ], what we have is the longitude given in Definition |2.2| which as we already stressed does not 
give a grading to the ring. You can readily check that for two units in different orthants |a + <x'\ < |a| + |a'|. 
For example, for D = 2 the Laurent polynomial L(zi,Z 2 ) = z 1 2 + z 2 2 + 1 is of that type, observing 
ziZ 2 L(zi,Z 2 ) = zi + zi + Z 1 Z 2 we deduce £(ziZ 2 L) = 2 < f(ziZ 2 ) + £(L) = 3. However, there are other poly¬ 
nomials that do satisfy this important property, for example L(zi, Z 2 ) = z\ 2 + z 2 + zi. Indeed, z“' z“ 2 L = 
z“ 1_2 Z 2 2 +z“ ,+2 Z 2 2 +z“ 1 Z 2 2+1 and i[ z*’z£ 2 L) = max(|ai — 2| + |a 2 |, |cxi +2| + |a 2 |, |ail + l «2 + 1|) = |a 1 | + |a 2 |+ 2 . 

Definition 3.7. A Laurent polynomial L is nice if 

£(LM) = £(L) + «(M), 

for every Laurent polynomial M G Clz^ 1 ,..., z^ 1 ]. 

Proposition 3.11. lf(V^,V^) C NP(L) with i (L) = m and a € Zd, then L is nice. 

Proof. Remember that 

£(z a L) = max (|a{+ ai|4-h |a[j + aol) < |a| + £(L), 

a'eNP(L)nz D 

and if a pair of opposed vertices are included in the Newton polytope, say V, m 1 , as max(|m + a.] |, | — m + 
ail) = |m| + |ai |, the bound is saturated 

£(z“L) = |a| + £(L). 


□ 


However, these Laurent polynomials are not the only nice Laurent polynomials, for example L 



+ 


Z 1 Z 2 1 + Z 1 Z 2 is also nice. Indeed, z“' z“ 2 L = z^'^zf 2 + z“ 1+1 Z 2 2 1 + z“ 1+1 Z 2 2+1 and t(z*'zf 2 L) = max(|ai — 
2| + |CX 2 I, |ai + 11 + \ot 2 — 11, |ai + 11 + |a 2 + 11) = |ai | + I 0 C 2 I + 2. Another example is L = izf ] Z 2 + 2ziZ 2 — 
iziz^ 1 + Izf^zf 2 + zi + zf 1 + 5, in this case we have t(z“'z“ 2 L) = max(|ai — 11 + |a 2 + 11, |ai +I| + |a 2 + 
11, |a] + 11 + |a 2 - 11, |ai I + |a 2 - 11) = |ai | + |a 2 | + 2. 

The following is taken from tropical geometry |65l 


Definition 3.8. (1) The Minkowski sum of the two sets A, B c M D is given by 


A + B = {a + b : a G A, b e B}. 


(2) For each w g M d zve can introduce an alternative order in lP and say a > w a' if and only if {a— a') xv > 
0 . 

(3) The initial form along the weight w G M D of a Laurent polynomial is 

in w L := f L a z , 

here the support is A w := {a G A : a' > w « Va r G A}. 

(4) The tropicalization of a Laurent polynomial also known as its tropical hypersurface is 

T{ L) := {w G M d : in w (L) is not a monomial}. 


Proposition 3.12. (1) For two Laurent polynomials K, L G C [z ± 1 ] we have the following formulae for their 

Newton polytopes and faces in terms of the Minkowski sum (the tropical product of polytopes) 

NP(KL) =NP(K) + NP(L), 


(3.11) 


F W (KL) =F W (K) + F W (L). 
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(2) The Newton polytope NP(z“L) and its faces are obtained from the Neioton polytope NP(L) via an automor¬ 
phism in the Abelian group Z D -a translation by a- 

NP(z“L) ={a} + NP(L), F w (z“L) = {a} + F W (L). 

(3) Newton poly topes of initial forms are the faces of the Newton poly topes 

NP(in w (L)) = F W (NP(L)). 


(4) The tropical hypersurface of a Laurent polynomial is 

T(L) = {w G R D : dimF W (L) > l}. 


TFie ortFiants of integers introduced in Definition 2.16 are extended to tFie real space as follows 


Definition 3.9. The real orthants in 


are 




i=1 



i G (a \ 3a), 
i G 3a, 

i G (Ca\ 3(Ca)), 
i G 3Ca, 


where R^ := R± \ {0}. For a = Z we define 

(M d ) Zd := (R_) D \ {0}. 

In this case we do not need to distinguisFi between R^ and its topological closure M,> 
consistency witFi previous developments we use a similar notation. 


M±, but for 


Proposition 3.13. Among the set of faces (F w (m) : dimF w (m) > 0} wgR D of the regular hyper-octahedron 
Con([m]) there are 2° facets {F <7 (m)} 0 . g2 * D/ with a facet per orthantF a [m) C (R D ) ff . Non facet faces of codimension 
p G {2,..., D} belong to the intersection p orthants. 


Definition 3.10. We introduce NP CT (L) := (NP(L)) n (R D ) ff . 


Proposition 3.14. A Laurent polynomial L is nice if and only zf NPo-(L) contains elements of longitude £(L )for each 
a G 2 Z °. 


Proof. See Appendix |C.6[ □ 

Proposition 3.15. If a Laurent polynomial L is nice then its tropical hypersurface T{ L) has a non trivial intersection 
with each orthant 

T(L)n(K D ) ff /l, Vug 2 Zd . 

Proposition 3.16. If a Laurent polynomial Lz's nice the polytopes corresponding to the translations by the generators 
oflP of its Newton polytope, {±e Q } + NP(L) have elements of longitude £(L) + 1 for each a G Zq. 

Proof. Let us suppose that a Laurent polynomial L is nice. Then, for each subset a C {1,..., D} we have a 
vector a a G (NP(L)) n (R D ) ff with |a ff | = £(L). Therefore, as e Q G (R D ) ff for a G Caand —e a G (R D ) a for 
a G a we have |a ff + e a | = |a a | + 1 for all a G Ca and |« ff — e Q | = |a Q | + 1 whenever a G a. Therefore, as 
this happens for each orthant, we conclude that ({±e a } + NP(L)) n {F w (Conv[f (L) + 1])} wgR D f 0. □ 

The previous two results are not characterizations of nicety but properties of Laurent polynomials that 
indeed could have non-nice Laurent polynomials. We now give a similar result to the last one that in fact 
is a characterization. For that aim we need to introduce some geometrical elements of the regular hyper¬ 
octahedron 


Definition 3.11. hi each orthant (R D )o- we consider the vectors 

Uc= Y_ e Q -^e Q . 

aeCa aGcj 
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Proposition 3.17. The facets of the regular hyper-octahedron Conv([m]) are given by F U(J (Conv([m])). 

Proposition 3.18. A Laurent polynomial Us nice if and only if{ u ff } + NP(L) has elements of longitude £(L) + D 
for each a e 2 Z °. 


Proof. =g For a nice Laurent polynomial in each orthant we have a vector such that £ (NP(L)) n (M D ) (7 
with | cxcrl = £(L). Therefore, as a a + u CT £ (M d ) ct so that |(Xo- + u a | = £(L) -F D. As this happens for each 
orthant, we deduce that ({u ff } -F NP(L)) n (F w (Conv[£(L) + D])} weRD f 0. 

4= We proceed by contradiction and suppose that there is no such point in a given orthant (R D ) a . Conse¬ 
quently, the set {u ff } + NP (r f L) contains only vectors of longitude less than £(L) + D. But this also holds for 
the rest of orthants component of the Newton polytope NIL' (L), cr a', as in that case, despite having vec¬ 
tors in the Newton polytope of the maximum longitude £(L), when we consider the set {u (I }+NI’ (I / (L) being 
the vector u ff and the vectors in NP a / (L) in different orthants, some of the components are summations but 
some are subtractions, and the bound £(L) + D is never reached. □ 

Proposition 3.19. Let L = Li L 2 be the product of two Laurent polynomials. Then, the Laurent polynomial L is nice 
if and only if Li, L 2 are nice. 

Proof. See Appendix |C.7| □ 

Corollary 3.1. The set of nice polynomials is a monoid under the polynomial multiplication. 

In the next picture we illustrate this situation, we have drawn Conv( [k]), for k = 0,..., 4, in dashed lines, 
and also the nice Newton polytopes NP(z^“ 2 +zi -Tzyzf], NP(iz 1 _1 Z2+2ziZ2—iziz^ 1 +2z 1 ~ 1 z 2 2 +zi -Fz^ 1 +9) 
and NI’(z 1 2 + z\ + zf] in different colors, red, brown and yellow, respectively 

a 2 



We now show the nice Newton polytopes and tropical curves in blue for basic degree 2 nice Laurent poly¬ 
nomials, any other one can be obtained from these by the corresponding convex hull of the unions of these 
basic Newton polytopes -observe that the two nice Newton polytopes of vertex type are not drawn here- 
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We show now some non basic nice Newton polytopes, and their tropical curves, built up by the convex 
hull of the union of the Newton polytopes of basic nice Laurent polynomials 




In the next picture we illustrate a case for D = 3, we draw two regular octahedra, the interior one corre¬ 
sponds to the convex hull of [ 1 ] and the exterior one is the octahedron corresponding to [ 2 ]; blue points are 
the multi-indices. Then, in red, we draw the Newton polytope of the nice Laurent polynomial, which is 
not a vertex pair polynomial, L = Z]Zt, -f 2z ] ] Zj, — Z 2 Z 3 — z 2 1 Z 3 + 3zJ 2 
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3.4. Laurent-Vandermonde matrices, poised sets and algebraic geometry. Here we follow | fl4| . The 
proofs can be easily translated from there. 


Definition 3.12. We consider the Laurent-Vandermonde type matrix 

v m ;= ( X W (pi)) ... )X M (pTi j) eC N Wm -AV ) 

made up of truncated vectors of multivariate units ^ k+Tn J(z) evaluated at the nodes. We also consider the following 
truncation S™ € (C rk ’ mXNk + m - 1 of the lower unitriangular factor S of the Gauss-Borel factorization of the moment 
matrix 


(3.12) 


cm_ 


( S M,[0] 

%+!],[< 


S M,[1] 

S [k+l],[1] 


%]| °[k],[k+1] 

S[k+l],[k] %k+1]| 


0 


[k],[k+m— 1] 


\ 


\S[k+m—1],[0] ^[k+m—1],[1] 

We have the following factorization 
(3.13) 
so that 


• • • 0[k+1],[k+ m —1] 

S[k+m— 1],[k+m—2] %k+m.—1]| / 


z* = srv k m . 


(3.14) KerV^ 1 C Kerl£\ ImlJJ 1 C ImS^ = C Tk -. 

The poisedness of j\4, m happens iff Ker X" 1 = {0} or equivalently iff dimlm I" 1 = r| Citll . 

Proposition 3.20. For a poised set A4,m the multivariate Laurent-Vandermonde V™ is a full column rank matrix; 
i.e., dim Im V™ = r k)m . 
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The study of the orthogonal complement of the rank; i.e, the linear subspace (imV ™) 1 C C Nklm 1 of 
vectors orthogonal to the image ImV™ where v G (ImV ^) 1 if = 0, will be useful in the study of 

poised sets. As Im V™ © (Im V™) 1 ^ = C Nk + m -' we have 

dim (Im V ™) ± + dim (Im V™) = N k+m _i. 

Proposition 3.21. The Laurent-Vandermonde matrix V™ has full column rank if and only if 

dim (Im V ™) -1 = N k+m _-, - r k)m = N k _i. 

Definition 3.13. Tor a Laurent polynomial K G C[zf ] ,... , z^ 1 ] the corresponding principal ideal is (K) := 
{z a K(z) : a G Z D }. Then, we introduce its truncations by taking its intersection with the Laurent polynomials 
of longitude less than k + m. We use the following notation for the truncated principal ideals 

(KW. := C{z“K(z) : a € Z D } [(z „ K(z||<k+m . 

It happens that the elements in the orthogonal complement of the rank of the Laurent-Vandermonde 
matrix are Laurent polynomials with zeroes at the nodes in the algebraic torus 

Proposition 3.22. As linear spaces the orthogonal complement of the Laurent-Vandermonde matrix (Im V™) 1 ^ and 
the space of polynomials of longitude less than k + m and zeroes at ATg m are isomorphic. 

Proof. The linear bijection is 

N k + m -)—1 

V = (w)Sr - 1 - 1 € (Imvr) ± <-> K(Z) = Y- v ^ z<Xi 

i=0 

where K(z) do have zeroes at A f™. Indeed, we observe that a vector v = (Vi)[i k 0 fm 1 1 G (Im V ™)" 1 can be 
identified with the Laurent polynomial K(z) = V[Z a ' L which cancels, as a consequence of vV™ = 0, at 
the nodes. □ 


Given this linear isomorphism, for any Laurent polynomial K with i (K) < k + m with zeroes at A f™ we 
write K G (Im V ™) ± . 

Proposition 3.23. Given a Laurent polynomial K G (ImVJ 1 ) 1 then 

(K) k+m _i C (lmVr) ± , 

or equivalently 

( K )k+m-1 ^ ImVf. 


Proposition 3.24. If a Laurent polynomial K is nice then 


and if K is not nice we have 


dim(K) k+£(K) _ 1 


N k _i, 


dim(K) k+£(K) _ 1 > N k _!. 


Proof. For a nice polynomial K we have t(z“K) = |a| + f(V) so that 

(K)k+i(K)—1 = C{z a K(z) : a G Z D }| a | <k 

and therefore 


dim(K) k+£ ( K )_i — |[k — 1]| + • • • + |[0]| — N k _i. 

However, if K is not nice then i(z a K) < |cc| + 1 ( K ) . Therefore, for all a G Z D with longitude less that |<x| < k 
we get £(z“K) < k + i (K), but there are also a G Z with |a| > k so that £(z a K) < k + £(K) and therefore the 
linear dimension of the truncated principal ideal is bigger than N k _i. □ 



34 


GERARDO ARIZNABARRETA AND MANUEL MANAS 


After this negative result, which shows that non nice Laurent polynomial perturbations of the measure 
are not well suited for the sample matrix approach to the Darboux transformation, we give a positive result. 
Following the ideas of |fl4| we find that nice Laurent polynomials are indeed nice with the sample matrix 
trick. 


Theorem 3.2. Given a nice Laurent polynomial L the Laurent-Vandermonde matrix V™ has full column rank if and 
only if 


(k)ic+«(L)-i — 


Proof. It is a consequence of the niceness of the Laurent polynomial L and Propositions 3.23 3.24 and 

MU □ 


From the spectral properties of the shift matrices we deduce 

Proposition 3.25. The row (L(Y)) W , a E Z D , is the longilex ordering of the entries in the corresponding 
polynomial z“L(z). 

Thus, in some way L(Y) encodes the same information as the principal ideal of (L) does. To make this 
observation formal we first consider 


Definition 3.14. We introduce the matrix (L(Y))^ k ’ m] E C Nk -' XTm - k given by 


( (L(Y)) [0]>m .. 

• (L(Y))[o],[ m _ic+i] 

0[o],[m-k] 


®[0],[k+m— 1 ] 

\ 

(L(Y))[ lc _ 2 ],[k] 


(L( Y) ) [ k —2], [k+m—2] 

0[k-2],[k+m-l] 


\(L(Y)) [k _ 1])[k] .. 

1 °[0],[k] 

°[0],[k+1] 

0[0],[k+m-1] ^ 


(k(Y))[ k _ 1 jj k+m _ 1 

/ 

0[k-m],[k] 

®[k—m],[k+1] 

• • • ^[k—m],[k+m—1] 




(h(Y)) [k—m+1],[k] 

®[k—m+1],[k+1] 

• • • ^[k—m.+1],[k+m—1] 

) 



(h (Y)) [k—m+2], [k] 

(h(Y)) [k—m+2],[k+1] 

0[k.—2],[k+m— 1] 




\ (L(Y)) [lc _ 1]>[k] 

(l-(Y))[ lc _ 1 ] ) [ k+1 ] 

••• (T-(Y))[ k _i] ) [ k+m _i ] y 





We collect this matrix and the truncation L(Y)^ in 

(L(Y))F := (L(Y) m , L(Y) [k ’ m] ) e C Nk -' xNk +™-'. 


Proposition 3.26. We have the following isomorphism 

K= Y- K a z“E (L) k+m _! <^(K ao ,...,K aNk+m _ i ) = (a 0 ,...,Q Nk _,)(L(Y))^ 

|a|<k+m 

between the truncated ideal (L) k+m _i and the orbit of C Nk -' under the linear morphism (LY))™.. Here zve have 
ordered the multi-indices a in L according to the longilex order, ao < • • • < «N k+m _! ■ 

Proof. Just recall that K is going to be a linear combination of the polynomials z“L(z), |ct| < k. Thus, the 
row vector by arranging its coefficients according to our order must in the rank of (LY))™. □ 


Now we have to show that full rankness of the Laurent-Vandermonde matrix is a sufficient condition 


Theorem 3.3. Let L E C[z ±] ] be a nice Laurent polynomial with £(L) = m. Then, the set of nodes C Z(L) is 
poised if and only if the Laurent-Vandermonde matrix V ™ has full column rank. 
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Proof. Let uvv s assume the contrary, then the sample matrix Z™ is singular, and we can find a nontrivial 
linear dependence among its rows (L™)i, i G {1,..., Tk )tn } of the form 

r k,m 

ci(4 n )i = 0, 

i=1 


for some nontrivial scalars {ci,..., c T|c m }. But, according to Lemma 3.13 (IjT)t = (S™)tV™, where (S™)i is 
the i-th row of S k and we can write 


. T k,i 




i=1 


so that 


T k, 

I 

i=l 


Ci (S^)t G (imV^ 1 ) 1 , 


and given the column full rankness of the Laurent-Vandermonde matrix, see Theorem |3.2j we can write 

(c l .---.c, k ,JSE'e(ImVr) ± , 

and following Proposition |3.26| we get 

(3-15) (ao, • • •, a N k _i) (L(T))™ = ( Cl ,..., c rki JS£\ 

for some non trivial set of c's. 

Now, we recall that S is lower unitriangular by blocks 



( S [k] I 0 

0 ... 

s = 

cm 

0 ••• 


V 


and that either L(Y) or L( J) are block banded matrices with only the first m block superdiagonals non zero 


L(T ) = ( (i-m) 1M (L(T)) 

l * * 


[k,m] 


0 


, {J) _( (L(I)) 1M (L(J)) 
l * * 


[k,m] 


0 


Now, we focus on the relation 

SL(Y) = L(J)S. 

which can be written as it follows 

(3.16) S [k] (L(Y))™ = (L(J)) [k,m] Sj^+ ((L(J)) [k] S [k] ,0). 

If we assume that ( 3.15| holds and left multiply ( |3.1 6| with the nonzero vector (a 0 ,..., a Nk , )(S k ') 1 we 
get 

(ao,...,a Nk _ 1 )(L(Y))^ = (ao,...,QN k _ 1 )(S [k1 )" 1 (L(J)) [k ’ m] S^ + ((a 0 ,..., a N]c _,)(S [k] y' (L(J)) [k] S [k] , 0). 

Thus, ( 3.15| holds if and only if 

(o 0 ,...,a Nk _ 1 )(S Ik] )- 1 (UJ)) M S [k] = 0, 

or equivalently if and only if det ^ (L( J)) ^^ = 0. Thus, recalling first Proposition 3.10 and second our inital 
assumption that the measures d p and its perturbation L d p do have MVOLPUT, the results follows. □ 


We say that a Laurent polynomial is irreducible if it can not be written as the product of non invertible 
Laurent polynomials —remember that the invertible Laurent polynomials are the units— . As the ring of 
Laurent polynomials C[z( 1 ,..., zff ] is a UFL^jany prime polynomial is irreducible. A Laurent polynomial 
L, different from zero or a unit, is prime if whenever L divides the product of two Laurent polynomials Li L 2 


'Up to units there is a unique factorization in terms of irreducible polynomials 
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then L divides L] or divides L 2 . Equivalently, L is a prime Laurent polynomial if and only if the principal 
ideal (L) is a nonzero prime idealfj 

Any Laurent polynomial L(z) = z“Q(z) can be written as the product of a unit z“ and a polynomial 
Q(z), and among the possible a there exists only one of minimum longitude^] To each irreducible Laurent 
polynomial L G Cfcsf 1 ,...,^ 1 ] it corresponds a unique irreducible polynomial Q G C[zi,... ,zq] and a 
unique unit az“, with L(z) = az“Q(z). Consequently, Z(L) = Z(Q) n (C*) D is the algebraic hypersurface 
Q(z) = 0 within the algebraic torus. 

Proposition 3.27. Any nice Laurent polynomial L could be expressed as the product of a unit u with different nice 
irreducible Laurent polynomials {Li,..., Ln} with multiplicities {di,..., dN} 

L — uL dl • • • L dN 

L — ULj l n 

in a unique form up to unities. 

Proof. Any Laurent polynomial can be factored in terms of its prime factors. But according to Proposition 
l3T9lif L is nice the factors must be nice as well. □ 

Proposition 3.28. Let L = uLi • • • Ln be a product of a unit u G Clz^ 1 ] with N different prime Laurent polynomials 
{Li,..., Ln), where none of them are units, and suppose that K G Cfc^ 1 ] is such that Z(L) C Z(K) then K G (L). 

Proof. Let us consider the product of a unit u G C[z ±1 | with N different Laurent irreducible polynomials 
{Li, • • • , Ln}, none of them units in Cfe^ 1 ], andlet{Qi,..., Qn) be the corresponding irreducible polynomi¬ 
als, none of them are units in C[z]; then, the related algebraic hypersurface is reducible Z(L) = U t N = 1 Z(Li) = 

(u t N =1 Z(Qi)^ n (C*) D . Now, given L = uLi • • • Ln let K € C[z{ 1 ,..., zf ] | be a Laurent polynomial such that 
Z(L) C Z(K), and write K = z“P where a has minimum longitude and P G C[zi,..., Z\f\ is a polynomial 
which is not a unit in C[z], Then, we know that ( Z(Qj j n (C*) D C Z(P), thus as Z(P), being the null 

set of a continuous function, is a closed set, we have that uJljZfQO = Z(Qi)) n (C*) D C Z(P). Thus, 

according to the Hilbert's Nullstellensatz P G -y/(Q), Q := nUi Qd being different prime polynomials this 
is a radical ideal and, consequently, P G (Q). Prom where we immediately conclude K G (L). □ 

Theorem 3.4. Let L = uLi • • • Ln be a product of a unit u g C[z ±] ] with N different prime Laurent polynomials 
{Li,..., Ln}, where none of them are units, and m = f (L). Then, the node set A/k,m C Z(L) C (C*) D is poised if the 
nodes dot not belong to any further algebraic hypersurface, different from Z(L), of a Laurent polynomial t with non 
trivial truncated ideal (t) k+m ^i 7 ^ { 0 }. 


Proof. Given that L is a nice Laurent polynomial we have that 

dirn(L) k+m _i = N k _i 

and as (L) k+m _i C (Im V,} 11 ) ^ we conclude that dim ( V™) ± > N k 1 . Prom Proposition 
are no more linear constraints derived from the inclusion of the node set in the algebraic hypersurface of 
L. Hence, if the nodes do not belong to any further algebraic hypersurface of a Laurent polynomial t such 
that (t) k+m ^i f {0} we have dim ( V ,} 1 ) 1 = N k _i and the set is poised. □ 


3.28 


we see that there 


3.5. Darboux transformations for the Lebesgue-Haar measure. As an example let us consider Darboux 
perturbations for the Lebesgue-Haar measure in the unit torus T° 


dp(0) 


1 

(27t) D 


d 0 i • • • d 0o 


9 An ideal in the ring of Laurent polynomials is a prime ideal if whenever the product of two Laurent polynomials Li L 2 belongs 
to the ideal we can ensure that at least one of the factors belongs to the given ideal. 

10 Take among all the mult-indices in the support of the Laurent polynomial those with the most negative component, say 
—At £ Z_, in each direction i € {1,..., D}, then can write L(z) = z~ A Q(z), and all the other units are of the form —A + ZP. 
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the MVOLPUT in this case are the units in the Laurent polynomial ring or one may say the multivariate 
Fourier basis 

= xik]- 

Now, the moment matrix is the identity matrix G = I, therefore the transformed or perturbed moment 
matrix will be 

TG = L(Y). 

Consequently, in order to have orthogonal polynomials, i.e. a Gaussian-Borel factorization, its block prin¬ 
cipal minors should be regular 

det((L(Y)) [k] ) /0, k G {1,2,...}. 

Which is achieved whenever the Laurent polynomial L is definite positive, i.e., 

det((L(Y)) [k] ) >0, k G (1,2,.. 

To illustrate this result let us explore a simple case, take for D = 2 the following nice definite positive 
Laurent polynomial L = zi + z ^ 1 + Z 2 + z^ 1 + 5, that preserves reality and definite positiveness. The original 
measure is the Lebesgue measure d0i d 02 in T 2 while the perturbed measure is (2cos(0i) + 2cos( 02 ) + 
5) d 0i d 02 - Now the longitude is f(L) = 1 and the Laurent polynomial can be written as follows 

L =z 1 ~ 1 zJ ] Q(zi,Z 2 ), Q :=zfz 2 T Ziz 2 + Z 2 + zi + 5 ziZ 2 

being Q G C[zi, Z 2 ] an irreducible polynomial (according to Maple). The corresponding real algebraic curve 
plot produced by Maple is 



Observe that for (zi, Z 2 ) G Z(L) we have 


1 


Z2 = -(-Zl-Zi — 5 ± F(zt )), 


-1 


1 


2 — 2 ( z i ^ T F(zi)) , 


where F := \J [z 1 + z 1 1 + 5) 2 — 4. 

We know that poised sets of zeroes of L do exist as long as we impose to these nodes not to belong to 
any other algebraic curve of a Laurent polynomial t having a non trivial truncated ideal (t)^. 

Let us give some details for the first non trivial case with k = 1; now we have ri j =4 and the set of 

2 ) t . The sample matrix is 


I = 


c (C*) 

where p, 

= (Pi,1 

fvT, I 

P2,1 

P3,l 

P4,l\ 

Pi,2 

P2,2 

P 3 ,2 

P 4 ,2 

Pi,2 

P2,2 

P3,2 

P4,2 

\P1,1 

P2,1 

P3.1 

P4,l/ 


G C 


4x4 
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Thus, if we request to the nodes to not belong to another algebraic curve as said before, we know that 
it is rank 4 and the set is poised. A similar conclusion is obtained directly, the determinant can be easily 
computed to be 


detl 


i _ 1 

1 “4 


= 5 


Pl.l 

-5 + F(pi,i) 
-5-F(pi,i) 

P 2,1 P 3 ,1 

-5 + F(p 2 ,i) —5 + F(p 3) i) 
-5-F(P2 ,i) -5-F(p 3 ,i) 

P4,1 

-5 + F(p 4) i) 
—5 — F(p 4) i) 

PU 

P2,1 

P3,1 

P4,l 

1 1 

1 1 



Pl,1 P2,1 

PiJ v 2 , 1 

P(Pl,l) P(P2,1 ' 

T 3‘ 

P- 

cff 

P- 



P3,1 P4,1 

) F(p 3)1 ) F(p 4>1 ) 




And we see, by direct substitution, that generically it does not cancel. 
For those sets of nodes the MVOLPUT of longitude 1 are 


4Ti](zi,z 2 ) 




/I 1 1 

0 0 0 

0 0\ 


1 


0 1 0 

1 0 1 

0 0 


zi + Zi ' + z 2 + Z-, ^ +5 

0 0 1 

0 1 0 

1 0 




o 

o 

0 0 1 

1 V 



( Pi,i 

P2,1 

1 ’ 1 

Pyi 

1 ’ 1 

P4~1 



P1,1P1,2 

P2,lP2,2 

’I ’ 

P 3 ,lP3,2 

P4,lP 4 ,2 

'1 ’ 


Pl,lPl,2 

P2,lP2,2 

P 3 ,lP3,2 

P4,lP4,2 


Pi,2 

P2,2 

F 3 ,2 

P 4 ,2 



Pi,2 

P2,2 

P 3 ,2 

P4,2 



Pi >1 Pi ,2 

P2,1P 2 J 

P3,1P 3 ] 

P4,lP 4)2 


P1,1P1,2 

P2,lP2,2 

P3,lP 3 ,2 

P4,lP4,2 


V Pi.i 

P2,1 

F* 3 ,i 

P4,1 

J 


(£) 


ZlZ^ 

Z]Z 2 


\ Z 1 

/ 





/ Pi ,1 
Pi, 2 

Pi,] 

P3,1 

Pa}\ 

-1 

f z A 

Pl,l 

P 3 ,2 

P 4 ,2 


z 2 ] 

Pi,2 

P2,2 

P3,2 

P4,2 


Z 2 

\Pl,1 

P2,1 

P 3 ,1 

P4,l/ 


\ Z 1 / 


Let us check that we are dealing with a vector of four longitude 1 Laurent polynomials. Observe 


(£A 


/I 1 1 0 0 0 0 0\ 


Z, ’z 2 


/Z, 1 (Z 3 1 +Z2 1 +Z 2 )\ 

0 10 10 10 0 




Z 2 1 K 1 + Z 2 1 + Zi) 

0 0 10 10 10 




z 2 (zf ] +Z 2 + Zl) 

0 

0 

0 

0 


Z1Z2 1 


V Zi (z 2 1 + Z 2 + Zi) / 


Zi Z2 

l 4 1 
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and 


1 0 0 0 0 0 \ 
10 10 0 
0 10 10 


/I 1 

0 1 0 
0 0 1 
\0 0 0 0 0 1 


1 1 / 


( ~- 2 


Pi,T P 2 ) 1 _ P 3 , 1 _ P 4 , 1 _ \ 

PyPli P 2 ,lP 2,2 P3,lP3,2 P4,lP4,2 

Pl,lPl,2 P 2 ) 1p2,2 P 3 iP3,2 P 4 jp4,2 
—2 —2 —2 _2 

Pi ,2 P 2,2 P3,2 P4,2 

Pl, 2 _ Vi, 2_ P 3 , 2 _ P4,2_ 

Pl,lPl] P 2 , 1 P 2 ,2 P3,lP3,2 P4,1P 4 ,2 

Pi ,1 Pi ,2 P 2 , 1 P 2,2 P3,1P3,2 P4,lP4,2 


V Pl.l 


P 2,1 


P3,1 


P4,1 / 


- 1\-1 


/-Pv(5 + Pl ,l) -P^{ (5 + P2,1) —Pi; 1 ! (5 + P3,1) 

~Pp2( 5 +Pl,2) -P^( 5 +P2,2) -P^ 2 (5 +P3,2) 

-Pl,2(5 + p^) -P2,2(5 + p^2) “P3,2(5+p^)-P4,2(5+P4^) 

Pi,i (5 + p7j ) -P2,1 (5 + Vi] ) -P3,1 (5 + Vi ]) -P4,i(5 + p^|) ) 


-p 4) ] (5 + p 4 ,i) \ 
-p 42 (5 + p 4 , 2 ) 


— —5 — 


ft )- 1 


/I 1 1 1 \ 
1111 
1111 
VI 1 1 1/ 


- 1\-1 


Then, we notice that from 


n i 1 n 


/! i i n 


n 1 i n 


n 1 1 i\ 

1111 

_ i _ 

1111 


1111 

/ _ i \ —i 1 

1111 

1111 

Li = —5 

1111 


1111 

P!) =5 

1111 

\i i i V 


\i i i V 


Vi i i V 


\i i i i/ 


we see that the marvelous event of the disappearing of the nodes occurs as we get 


n 


i 

0 1 

0 0 


0 0 0 0 o\ 


1 0 

0 1 


Vo 0 0 0 0 


0 0 

1 0 
1 1 


< h. 

Pi ,1 Pi,2 P2,lP2,2 P 3 ,lP3,2 P 4 ,lP 4 ,2 

'1 ' '1 ' '1 ' '1 ’ 


2 

-,1 

C -1 


Plj 
Pi ,2 


V pf,i 


P 2,1 

-C-i 


P3,l 

C -1 


P 2,2 

P 2,2 


pJ 

P 3 ,2 


pii 


P3,l 


P4,1 

-C-i 


-i. 


Pi,1 Pi,2 P 2 iP2,2 P 3 iP3,2 P 4 iP4,2 


P 4 j 

P4,2 


Pl,lPl, 2 P 2 , 1 P 2 , 2 P3,1P 3)2 P4,iP 4i2 
Pi ,1 Pi ,2 P 2 , 1 P 2,2 P3,1P3,2 P4,lP4,2 


P4,l 


/ 


A\] P 2 ,i 
Pi ,2 Pli 

Pi ,2 P 2,2 

Vpi ,1 P 2,1 


P3,l P4,l\ 
P 3 ] P4,2 
P3,2 P4,2 

P3,1 P4,l/ 


n i 1 a 


= - 5+ 5 


1 1 1 / 
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Moreover, we get 


/I 

0 

0 

Vo 


°\ 

0 

0 

V 


( Pij 
PpiPlj 
P1,1P 1)2 
Pij 
Pl,2_ 
Pl.lPlj 
Pi ,1 Pi ,2 

) Pl.l 


P2,1_ 

P2,lP2,2 

P2,]P2,2 

-2 

P2,2 

2 

P2,2_ 

P2,1P 2) 2 

P2,lP2,2 

pli 


Ps,1 

P3,lP3,2 

P3,]P3,2 

P3,2 

P3,2_ 

P3,1P 3 ,2 

P3,1P3,2 

P3,1 


P 4 j_ 

P4,lP 4 ,2 

P 4 ,1P4,2 

—2 

P 4 ,2 

7 

P 4 ,2_ 

P4,1P 4) 2 

P 4 ,lP4,2 

P 4 ,1 


K\ 


(*!) 


U-l 


.-1 


Z2 

\zi ) 




n 

1 

1 




( z A 


^ 5z 1 1 + 1 — ^ 


^(zi + 5) - ^ 


5 1 

1 

1 

1 

1 



Z , 1 


5z 2 1 + 1 - \ 


z 2 1 (z 2 + 5) - 5 


5 

1 

1 

1 

1 



Z 2 


5z 2 + 1 — | 


z 2 (z 2 ] +5) - 5 



V 

1 

1 

V 


\Zl / 


v 5z i + 1 - 5 ) 


W 2 ^ + 5 )-5/ 

Now, we come back and explicitly evaluate 4>[i] to get 






4>[1] = j 


/z, '{z ] 1 +z 2 1 + z 2 )\ 


/s 


1 + Z 7 1 + Zl _ 


Z2( Z-[ +Z 2 +Zl) 

V Z] (z 2 1 + Z 2 + Zl) J 


+ 


z 2 {z 2 


1 +S)-5 


z 7 (z 2 + 5) - 


V z ! (z, 1 + 5) 


y 


/V - i\ 



ZJ 1 - J 



Z 2 5 

/. 


V z l - 5 / 


and we see that another marvelous event occurs, the Laurent polynomial in the numerator factors out 
having as one of the factors the Laurent polynomial in the denominator which can be cleared out. Finally, 
we get the longitude one Laurent polynomials, orthogonal to 1 according to the measure L d 0i d 0 2 . 

For k = 2 now we have r 2 ,i = 8 and the set of nodes is A/ 2,1 = {pi, p 2 > P 3 > P 4 > P 5 , P 6 > P 7 > Ps) C Z(Q) c 
(C *) 2 where p : = (pi 


A = 


Pij 

Pi,I Pij 

p 1,1 pi >2 

Pij 
Pl,2_ 
Pi >1 Pi ,2 
Pl,lPl,2 
Pl,1 


L, 2 ) T - The 

sample matrix is 





—2 

-2 

-2 

-2 

-2 

-2 

-2 

P2,1 

P3,1 

P 4 ,1 

1 ’ 1 

P 5 ,1 

Pe,i 

1 1 1 

P7,1 

1 ’ 1 

Ps,i 

1 ’ 1 

P2,]P2,2 

P 3> 1 P 3,2 

P 4 ,lP 4 ,2 

P5,1p 5 ,2 

Pe,iP6,2 

P 7 ,lP 7 ,2 

P8,lP 8 ,2 

'-I ' 

P2,lP2,2 

’ -j 

P3,lP3,2 

P 4 lP4,2 

' -J 

P 5 ,lP5,2 

P 6 ,1p6,2 

' -J 

P 7 iP 7 ,2 

P8lP8,2 

P2,2 

P3,2 

P 4 ,2 

P5,2 

Pe,2 

P 7 ,2 

P8,2 

P2,2 

P3,2 

2 

P 4 ,2 

' I 

P5,2 

' -I 

2 

Pe,2 

' -I 

P 7 ,2 

pi ,2 

P2,1P 2 ,2 

P3,1P 3 ,2 

P4,1P 4 ,2 

P5,1P 5 ,2 

P6,1P 6 ,2 

P7,1P^2 

P8,1P 8 ,2 

P2,lP2,2 

P3,lP3,2 

P4,lP4,2 

P5,lP5,2 

P6,lP6,2 

P 7 ,lP 7 ,2 

P8,lP8,2 

P2,1 

P3,1 

P 4 ,1 

Ps,i 

P6,1 

P 7 ,1 

P8,1 


G C 


8x8 


and corresponding La uren t-Vandermonde matrix was explicitly above. Let us illustrate the argument in 
the proof of Theorem 


3.4 


We have that all nodes are solutions of L = 0 that is 5 = — fzi 


1 + z 2 ^ + z 2 + 


zi). Moreover, if we consider the ideal of L and write down the corresponding relevant equations at this 
moment, i.e., z-\L[z-\,z 2 ) = 0, z 2 L[z],z 2 ) = 0, zj'L[z-\,z 2 ) = 0 and z^L( zi,z 2 ) = 0 we get the following 
relations 

1 + 5zJ ] = — (zf 2 + zj ] z 2 + z^zj 11 ), 

1 +5z 2 ] =— (z-\z 2 ] +zj ] z 2 ] 1 ~~ 2 '' 

-1~- 1 


< 


z^L(zi,z 2 ) = 0, 
zf 1 L(zi,Z 2 ) = 0, 

z 2 L(zi,z 2 ) = 0, 
ziL(zi,z 2 ) = 0, 


< 


'- 2 T z 2 j, 


1 +5z 2 = — (Z 1 Z 2 + Z 1 ' Z 2 + Z 2 j, 
1 +5zi = —(z 2 + z-\z 2 + zi z^ 1 ). 


From these observations we see that the matrix V] is structured in 3 bands, being the last one our X 2 , that 
as the nodes belong to the algebraic hypersurface of L, give us 5 linear relations among the rows. That is, 
the rank is at much 8 . From previous results we know that if the set of nodes is poised then the rank is 8 . 
But observe that there are no linear dependence relations among the rows in Z 2 , all the 5 relations involve 
rows in the previous bands, not in the last one L\. 
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3.6. Perturbed Christoffel-Darboux kernels. Connection formulae. We discuss here how the Christoffel- 
Darboux kernels behave under Darboux transformations, expressing the perturbed kernels in terms of the 
non-perturbed ones 


Definition 3.15. We introduce the following tzvo truncations of the resolvent matrix. The first one is an upper block 
triangular matrix 



(cv 


tu [l],[H-1] • 

“>[l],[l+m-1] 

\ 

tu [l ’ m] := 


0 

“>[1+1],[1+1] • 

“>[1+1],[l+m] 



V 

0 

0 

“> [l+m— 1 ], [l+m- 

-1 ]/ 


ivhile the second one is a loiver block triangular matrix 


( “>[l],[l+m] 

0 

.. 0 

\ 

“>[l+1],[l+m+1] 

“>[l+1],[l+m+2] 

0 


\“>[l+m-1],[l+m] 

“>[l+m—1],[l+m+1] 


a, [l+m-1],[l+2m—11/ 


and also the following block diagonal matrix 


H [l ’ ml := diag/H^,... ,H [l+m _i]) 


Notice that when cu l,m and cu[ l m ] are glued together we get a m block wide horizontal slice of the 
resolvent uo. We will also use cu L ^ 11 m which denotes the truncation of the resolvent cu built up with the 
first l block rows and the first l + m block columns. 


Theorem 3.5 (Connection formulae). The following formulae relating Christoffel-Darboux kernels hold true 


K (l+m) (z!,z 2 ) = L(z 2 )TK' l+Ta) ( z i, z 2 ) - 

/ T$[ l ](z 1 ) \ 

T$[i+i]( z i) 

t 

(7H [l+m ) m] ) -1 a , [iim] 

( 4* [l+m] ( z 2) ^ 

4>[l+m+1]( z 2) 


\T$[i+m—1] ( Z 1)/ 


\4>[l+ 2 m—1 ]( z 2)/ 



f T$[i](zi) \ 

t 

/ <t>[l] ( z 2) \ 

L(Z 2 )TK (1) (Z!,Z 2 ) + 

T$[i+i]( z i) 

(TH [l > m] ) -1 u) [l,m] 

4*[l+1] ( z 2) 


\T$[i+m—1] ( Z 1)/ 


W +m -p(z 2 )J 


= L(z 1 _1 )TK- l+m * (zi, z 2 ) 


( $[l+m]( z l) \ 

$[l+m+1]( z l) 




/ T4) [:] (Z 2 ) \ 

T4)[ l+1 ](Z2) 


\$[l+2m—1] ( Z 1)/ 


\T45[l+m-1]( z 2)/ 



/ ^$[i]( z i \ 

t 

( ^ m (z 2 ) \ 

Uzf ] )TK^{z u z 2 ) + 

$[i+i]( z i) 


T4>[l+1](Z2) 


\$[l+m—1] ( Z 1)/ 


\T4>[l+m-1]( z 2)/ 


Proof. First we consider the expression 

[l+m] 




(K/l+FT 1 


4) [l ’ ml (z 2 ) = (T$ [l ’ m] (z 1 ))t 


f(TH- 


\ [l+m] 


jO) 
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which of course is an identity due to the relations between A and to. We now let the operators act to their 
right or to their left and obtain 

(T$(zi)) t (rPi t hr i y 1 4>(z 2 ) = (T$(z 1 )) t ((Thr 1 )^) 1 1 «|>(z 2 ) . 

The term between brackets on the LHS is 

($^(z 1 )) t (H _1 ) [l+m] 

While the term between the brackets on the RHS is more complicated. Observe that 



To conclude let us say that the third and fourth equalities follow similarly from the equality 

($(z 1 )) t (TT 1 M) [l+mJ T4 >(z 2 ) = ($(zi)) f T4 >(z 2 ). 


4. INTEGRABLE SYSTEMS OF TODA AND KADOMTSEV-PETVIASVILII TYPE 

The idea here is to introduce deformations of the moment matrix or, equivalently, the measure. These 
deformations will depend on some continuous and discrete parameters. For those cases where the initial 
moment matrix is Hermitian (and therefore the measure is real) and/or definite positive we give conditions 
under which deformations preserve these properties. 


4.1. Discrete Toda type flows. Let us consider the composition of two Darboux transformations, for that 
aim we collect what we know of each single Darboux transformation as follows. Given two Laurent poly¬ 
nomials Li (z) and L 2 (z) we consider the corresponding Darboux transformed measures 

Tidp(0) = L i (e i9 )dq(0), 1 = 1,2, 

for which we have the resolvent and adjoint resolvent matrices 

oh = (TiS)Li(T)S _1 , d>t = (TiSMYjS- 1 , 

M i = ST i S-\ At = $Ti§-\ 


connected by 


hkbJTiHr 1 = Mi, 


(TiH) 1 uhH = 
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and such that the associated MVOLPUT satisfy 

oh®(z) = Li(z)Ti<D(z), d>P(z) = )Ti$(z), 

MiTiO(z) = ®(z), MiTid?(z) = 6(z). 

The composed Darboux transformation corresponds to the perturbed measure 

TtTi d p(0) = T^ d p(0) = Li (e i0 )L 2 (e ie ) d p(0). 

Observe that if we want to apply the sample matrix trick to find the resolvent we must require Li and L 2 
to be nice, and as was proven in Proposition [3T9] the product 0 L 2 is nice again. Therefore, the iteration of 
two nice Darboux transformations is a nice Darboux transformation. The interplay of these two Darboux 
transformations can be described by 


Proposition 4.1. The Jacobi matrices J Q , a = 1,..., D, fulfill the discrete Lax type equations 

(TiJa)tUx = CUiJ Q , M|T|J a = JqUH, 

for i = 1,2 and a = 1,..., D. Moreover, the following Zakharov-Shabat compatibility equation is satisfied 
(4.1) (Ticu 2 )uq = (T 2 cu-| )cu 2 , M] T| M 2 = jVl 2 T 2 ]VLi. 

The hatted versions of these equations also hold. 

Proof To prove the first set of equations just replace J Q = ST U S 1 and the expressions for cut and given 
in ( |4.1[ ). As for the second one, we perform the very same replacement to get 

M 2 - SL 2 (T)Li (rX(TiT 2 S)-\ 

from where the conclusion follows immediately from the symmetry in the interchange of subscripts. □ 


Therefore, T| and T 2 could be understood as two compatible discrete flows introducing in this way two 
directions for the associated discrete shifts. We can generate a lattice of discrete flows by considering d, in 
principle not related with D, possible directions or family of transformations or perturbations. For each 
i E (1,..., d} we need an infinite family of multivariate Laurent polynomials {L lJU (z)} ni£ z and consider a 
measure that depends on the multi-index m = (mi,..., m d ) e Tf\ as follows 

d Pm+ei(0) = Ti d p m (0) := Li >mi +i (e ) d p m (0). 


A simplified version of this lattice of measures can be gotten with 

d 

(4.2) dp m (0):=n(Li(e ie )) mi dp(0). 


i=1 


Notice that the introduction of d different directions is, in principle, arbitrary, and that we coud even 
collapse all to one direction d = 1 having the very same measure. This is connected with the bijection 
existing between and Z + . From the viewpoint of the Gauss-Borel factorization the moment matrix for 
Mis 


G(m) = n (Um) mi G = G Yl (Limp. 

i=l i=l 

We now consider discrete flows associated with Darboux transformations generated with degree one 
Laurent polynomials. Our discrete deformations needs . 


Definition 4.1. Let n Q = (n Q _i,..., u Q) ^d> n. Q) D,..., n Q) i ) T E C 2D be 2D linearly independent vectors and let 
q Q , a E {1,... ,2D} and m E Z 2D be complex numbers and integer multi-indices, respectively. Then, the evolved 
measure is 

2D 

L(m,z) := P[(L n Jz) - q Q ) ma , 

a=1 


d p m (z) :=L(m,z)dp(z), 
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the vacuum wave matrix is 

2D 

W 0 (m) :=[]^a-qa) m ‘ 

a=1 


and the deformed moment matrix is G(m) = Wo (m)G = GWo(m). 
Proposition 4.2. For degree one discrete flows reality is ensured if 

q a £ Ti a = n Q , 

This choice also ensures the nicety of the Laurent polynomial L(m, z). 
Equivalently, 


Va E Zd. 


UQ - ftQ - (Suuj' 


u Q E C 


D 


We see that for discrete transformations we have an one step real reduction, however we could also 
consider as in IfTlI two step reductions, in that case we request for 


^2DTIq = TId+1-q = j > 


qD+i-a = q a> 


m.D+1—a — TUaj 


which naturally appears when we ask for (L na (z) — q a ) mi • (L nD+la (z) — qD+i-a) mD+1_a to be real; involving 
therefore two products of Wo. Now the action of only one T a does not preserve the Hermitian character, 
but the composition of translations T a • To + i_ a does. 

The definite positiveness for the case considered in Proposition [42] is achieved if 

D 

2^ |u a , b |cos(e b + argu Q)b ) ^ q Q , a E {I,...,2D}, 

b=l 


which is fulfilled whenever 


D 

qa T 2 ^_ |LL Qj bl ^ 0. 

b=l 


Notice that 

(4.3) T a G = [yL a — q a )G = G()f a — q a ). 

with corresponding degree one resolvents given by 

o>a := (T Q S)(Tla - qa)S~ 1 > d)t := (§- 1 ) t ( T t a - q Q )(T Q §)t. 


For these adjoint resolvents we have 
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Proposition 4.3. The follozving expressions for the degree one resolvents hold 
a) a = (T Q H)(§T Q §- 1 ) t H- 1 

/(T a H[ 0] )Hjg (Tt a ) [0])[ i] 0 0 

0 (TaH^H- 1 (TtaJm.Ki 0 

0 0 (T a H [2 ] )H [2 j I (?ia)[2],[3] 

0 0 0 (TaH^H- 1 

V 

=H- 1 S(T a S)- 1 T a H 

/HfojTaHio] 0 0 0 

Wa)[i],[0] H^T Q H n] 0 0 

o (TU[2][,1] H [2] lT a H [2] 0 

0 0 (TU[3],[2] H [3] 1T Q H [3] 

V 




Proof It is a direct consequence of the first relations of T Q G = [yi a — q Q )G and the second of T Q G = 
G[yi a — q Q ). □ 


It also implies the following relations 

(T a H [0] )H^ = -(j* Q )[o],[i]0[i] - q a , 

(TqHmJH^ = (T a (3[ k ])('pf a )[ k _ 1 ] i [ k ] — (pf a )[k],[k+i]P[k+i] — q Q , 

(T a H[ k ] )(7f a )[ k ],[ k +i]H[ k+ i] = —A a |3[ k+1 ], 

H [0] (T a H)[o] = — P[i](l^a)[1],[0] qa> 

= Wa)[k],[k-1]Ta3[ k j — P[k+1] Wa) [k+1],[k] — qa> 
hd[k] (ff-a)[k],[k— 1]T a H[^_ij = A a (3[ k ], 
from where we derive 

Corollary 4.1. Discrete Toda type equations are fulfilled 

A b ((T Q H [k] )H-;) = [yL a ) Mjk+u'dlk+i] Wb)[k+i],[k]TbH [k] 1 - (T a H[ k] )(]l b )[ k ] )[k _-,](T a T b H [lc 1 _ 1] )(-pf a ) [k _i] ! [ k ], 

or equivalently 

^ b ( H [ k ]T a H [k]) = (TbH [k ] ) (ft b )[ k ] t [ k+ !]H [k+ 1J (tf Q ) [ k+ !] ( [ k ] - ("pta) [ k ],[k-l] (T a TbH ) [yi h ) [lc— 1 ] , [ k ] T a H [ k ] , 

for any couple of indexes a, b 

4.2. Miwa type expressions for the MVOLPUT. We introduce the matrices 

a a,[k] : = (T Q H[ k ])Hj k |, ^Q,[k] := ^[k] ^ahl[ k ]. 

In terms of which the nonzero degree one resolvent block entries are 

(<^a)[k],[k] = a a,[k]> (^a)[k],[k] = &a,[k] 

on the main diagonal, and (pf a )+ for cu Q and (ft a )+ for <£> a , on the first superdiagonal. The adjoint resolvent 
matrices are defined by 

At := T q H _1 cu q H = ($(T a Sr')\ M q := Hd^H" 1 = S(T a S)-\ 

and consequently can be written as follows 
M a = I + H(pf Q )_T Q TT\ 


At =I+(T Q Hr 1 (pt Q ) + H. 
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Hence, these block unitriangular matrices have only two block diagonals different from zero, the main 
diagonal which equals the identity and the first sub-diagonal with 

( M a)[k+1],[k] = Pa,[k+1]> (^a)[k+1],[k] = Pa,[k+1]> 


where 


f>Ul] ( TaH [k] )Wa)[k],[k+1] H [k+1]> 

Pa,[k+1] := Hnc+llWailk+ll.IklTaH^]. 

As we know, we have for the transformed and non-transformed MVOLPUT the following connection 
formulae 


CU a <X>(z) = (L na (z) - qJTa^z), 
A Q T a d?(z) = d(z), 

Now we need two technical lemmata 


<&a$(z) = (Lft a (z) - q a )T a d?(z), 
M q T q cD(z) = O(z). 


Lemma 4.1. The inequality 2D | [k] | ^ | [k + 1 ] | holds true for all k = 0,1,_ 

Proof. Given a multi-index a = (oq,..., <xd ) £ [k] such that its i-th component satisfies oq ^ 0, then 
a ± ei E [k + 1], Therefore, for each a E [k] and each of its D components we can get one (oq f 0) or 
two (at = 0) multi-indices in [k + 1] out of the initial a E [k]. Hence, it is clear that given any multi-index 
a E [k] we can get n multi-indices in [k + 1] with D ^ n ^ 2D. Now, we observe that any a' E [k + 1] can 
be obtained by this procedure form at least one multi-index in [k]. Then u|[k]| ^ |[k + l]|and consequently 

2D|[k]|^|[k+l]|. □ 


Lemma 4.2. The matrix 


(( T 1 1 )[k],[k-(-1]\ 

OHz^fcUk+l] 


(X[l](' Y '))[ k ] ! [ lc+1 ] 


(^D 1 )[k],[k+1] 
("^D) [k],[k+1] 


G M 2D|[k]|x|[k+1]| 


V fY'l)[k],[k+1] / 


has full column rank for all k = 0,1 .... 

Proof. See Appendix |C.8| 

Definition 4.2. For the degree one discrete flows zve consider 
(1) The matrix 


□ 


/n? 


T ’ 


N := 


E C 


2Dx2D 


V n 


-2D/ 


which zve assume to be invertible (as a consecuence N = NT 2 D is also be invertible). 
(2) For any vector z E C D the associated vector 


q := Nx[ij(z) E C 


2D 
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(3) We consider two rectangular matrices 


( (*i) 




[TH] k := 


[k],[k+1] 


-<2D|[k]|x|[k+1]| 


\ (l^2D) [k], [k+I ] / 

/TiH M 

\T2DH[ k ] / 


c 2D|[k]|x|[k]|_ 


We are ready for the following result expressing the MVOLPUT as products of the quasi-tau matrices H 
and its discrete flows 

Theorem 4.1. The MVOLPUT can he expressed in terms of quasi-tau matrices H and its degree one discrete flows 
as follows 

4>[k] (2) =(-1 [ TH ]k-i^ k _ 2 [TH] k-2 Hj^l_ 2] • • • , 

$[k] ( z_1 ) =(-l ) k $k-i [TH] k _!Hj^ # J_ 2 [TH] k _ 2 H^ 2] • • • $ [THJoH^. 

Proof. Since N is invertible the matrix 1)4 k still has full column rank. Hence, the correlation matrix 1)4 k f>4 k G 
Cl[k +1 ]|x|[k+lll i s invertible and the left inverse of 1)4 k is 

wt ■■= r 1 ^. 

All the previous said, and as long as the two previously stated conditions hold we have 

a [k]4>[k]( z ) + ( n a • X[l] ('^'))[k],[k+l]^[k+1] ( z ) = 0) 

whenever the transformation parameters {n Q , q a }^ are as described in Definition |42] and therefore 

^k^l43[k](2) = -[TH^H^cJy^z) 

But, since h k -i has full column rank |[k]| with left inverse 1)4^' , we get 

4>[k](z) = -h>J c _ 1) [TH] k _ 1 H^l_ 1] (l) [k _ 1] (z). 

It is easy to see that a similar result can be obtained for $ when evaluated z 1 . Iteration leads to the desired 


result. 


□ 


4.3. A Toda type integrable hierarchy. For the continuous deformation we introduce a covector t := 
(t[ 0 ], t[i],...) where t[ k ] = (t ( k ),..., t (k) ) are complex covectors as well. 


1 |[k]| 

Definition 4.3. The vacuum wave matrix is 

W 0 (t) := e t[z \ 

Then 


t(z) := t a T a - 

a€Z° 


Proposition 4.4. The matrix G(t) = Wo (t)G = GWo(t) is the moment matrix of the deformed measure 

d p t = e 1 ^ d p, 

and the Fourier series of the evolved measure is 

At = e t(z) (1. 


Definition 4.4. We say that t is an admissible set of tunes if G (t) admits a block Cholesky factorization as in 
Proposition \2.4\ 

Hereon we assume that we are dealing with admissible times. 

For the continuous flows we need exp(t(z)) to be real which is achieved when t(z) € R, this last require¬ 
ment holds for t such that t = tq; then, we request t[ k ] = t[k]£|[k]| or, equivalently. 
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Proposition 4.5. The continuous flows preserve the reality and definite positiveness whenever = t„. The real 
deformed measure in T D has the following form 


7 V“oo v— I [k] I /2 , 

z 2_k=o2-j-’ z 


, I. (ic) cosfa! kl -0} 

a > 


dM9)=e “J d p(0). 

From the Gauss-Borel factorization of the perturbed measure 

(4.4) W 0 (t)G = GWo(t) = G(t) = (S(t))“ 1 H(t)((§(t))- 1 ) t 
we get t-dependent MVOLPUT. 

Definition 4.5. The wave matrices are defined as follows 
W, :=S(t)W 0 (t), 

Wj := (S(t)) 1 H(t), 

Proposition 4.6. We have 

(4.5) G = Wf’Yfyz = W+ (^f 1 
Proof. Use the Gauss-Borel factorization of G(t) = Wq (t)G = GWo(t). 


W| := Wo(t)(S(t)) T , 

Vt' 2 :=H(t)((§(t))- 1 ) t . 


□ 


The splitting as a direct sum of the linear semi-infinite matrices in strictly block lower triangular matrices 
and upper block triangular matrices is denoted by M = M< + M>, and the associated splitting in strictly 
upper triangular and lower triangular parts is denoted by M = M< + M>. 


Lemma 4.3. We have the Gelfand-Dickey expressions 

+st“s _1 = (sr a s- ] 

Ot(x V 

4 ?-§^+§ T “$- 1 = / § Y a § _1 

3t_« V /> 


9S S” 1 = -( ST a S _1 


at 


=-f§T“§- r 

3t_ a V /< 


□ 


Proof See Appendix |C.9| 

Observe that for Hermitian cases the first and second lines in the previous lemma coincide —due to the 
equality between hatted and non hatted terms and t a = t a . 

Definition 4.6. We introduce the Zakharov-Shabat matrices 

B«:=(J«(t))>, := (T«(t))>. 

The next important result collects the essential elements of the associated integrable theory 

Proposition 4.7. (1) The wave matrices solve the linear systems 


(4.6) 


(4.7) 


3Wi 

3t« 

3\^ 2 

3t„ 


=B«W 1 , 

=Ba^ 2 , 


3Wi 

3t_« 

3W 2 

3t_„ 


=B«W! 


=B«W 2 


(2) The Jacobi matrices are Lax matrices ; i.e, they fulfill the following Lax equations 


3Jq 

3t„z 


— [B«', Jot]> 


3 Jc 


3t_ 




(3) The Zakharov-Shabat (or zero-curvature) equations are satisfied 


3B, 


3B r 


3t a / 3t a 

Proof. See Appendix |C.10| 


+ [B«, B a /] —0, 


3^ a 36 a / 


3t—a' 3t- a 


T ^a'] —3- 


□ 
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Notice that only the equations on the left or the right would be need in the Hermitian case. 
We now introduce the Baker functions and their adjoints 

Definition 4.7. The Baker and adjoint Baker wave functions are 


¥1 :=W lX , 
^2 := V^X, 


t 


% := (wr 

%= (W)'x, 


0 ) x, 


:= X^f, 
^ := X f WJ, 


(^) f :=xW> 

:=X^. 


Proposition 4.8. The Baker functions and adjoint Baker functions in Definition ^. 7\ can be expressed in terms of the 
MVOLPUT and the Fourier transform of the measure asfollozvs 

=e t(z) 0(z,t), z E (C*) D , ¥t = (27t) D a(z- 1 )((H(t))^ 1 ) t a(z,t), 

#1 = a(z,t), z E (C*) D , ^ = (27t) D jl(z)(H(t))^ 1 0(z,t), 

= (27t) D p:(z)e t(z) 0(z,t), ze®^, ^ = ((H(t)) -i y6(z,t), 

¥ 2 = (27t) D A(z- 1 )e^ 1 ^a(z,t), z^'e®, M^ = (H(t)) _1 ©(z,t), 

Proof. We use Definitions |2.10} [4.5| and [47| together with Proposition |2.10| 

Proposition 4.9. The Baker functions satisfy the following linear differential equations 

9¥i 


z 1 E V^, 
z E T>^ 
ZE (C*) D , 
ZE (C*) D - 


□ 


3t« 

9l>2 

9t„ 


=B«¥i, 

=bA 


9^1 
3t_ a 
9HL 
9t ~ 


=$0^2 


while the adjoint Baker functions satisfy 


9^>* 

3t« 

9Md 




9¥ 




9t-« 

9^: 


L =-(B_ a ) t ^, 


- =-(B_ a ) t ^. 


9t« v 7 9t_« 

Proof. Taking derivatives on their definitions and using Lemma [473] the result is straightforward. □ 

We will be specially interested in the tpj part of the covector of the continuous deformation parameters 

t[1] = (t_i,...,t_ D ,t D ,...,ti). 

Proposition 4.10. The following differential equations are satisfied by the square matrices H[ k ] and the rectangular 
matrices |3[ k ]. 


9H 


3p[k] 1_J f'y ^ U“1 

-0^“ = — H-[k] 1 ' aJ[k],[k-l] ^[k-lj 

M u _i 


3p[k] 

9t_ a 


**- Ht M (T Q ) H ,Dc-il (H 


-l 

[k-l] 


H [k] = P[k](^a)[ k -1],[lc] - (Y a )pc],[ k -|]P[ k +1] 


9H„ 


Hfoi 1 =-(Ta)[0][1]|3[1] 


9 t [k] l“Lltj • u>LK.-IJ,LK.J V • u n K.j,LK.-IJI“LK.-|-IJ gj. 


Proof. These relations are just identifications if we use in Lemma [43] the explicit formulae (2.141. 

Taking second partial derivatives in the previous expressions leads to the following 
Theorem 4.2. The H[ k ] matrices are subject to the following Toda lattice type equations 

9 ( 3H [k] 


□ 


(48) hr 

9tb 


9t a 


H [k J ) = (T Q ) [k])[k+1] H [k] (T b ) [k+1])[k] H [k j - H [k] (T b ) [k])[k _ 1] H [k 1 __ 1] (T a ) [k _ 1]i[k] 
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which can he rewritten for the rectangular matrices (3fy as 

9 2 |3ry 3 3p|kl 3(3r k i 

9tn g tb = gjr (P[k]( T b)[k-i],[k]P[k]) - -g^-|3[k-i]('Y'b)[k-2],[k-i] - ( T b)[k],[k+i](3[k+i]-^ tQ - 

andfor the as 

9t_ a ai_ b = WT a (^M( T b)[k-1],[k]£[k]) ~ g^-^[k-1](^b)[k-2],[k-1] - ('Y'bJM.Dt+ii^Dc+H g^- 

4.4. Bilinear equations. We begin with the following observation 

Proposition 4.11. Wave matrices evaluated at different admissible times t, t 7 are related as it follows 

Wi (t) (W, (t')r 1 = V^ 2 (t) (V^ 2 (t 7 )) -1 , (t) (ft, (t 7 )) -1 = W 2 (t) (W 2 (t 7 ))- 1 , w, (t) (w 2 (t , )) t = ft 2 [t) (ft, (t 7 )) 1 . 

Proof. Just remember |4.5| wich is valid when evaluated at any time, therefore 

G = (W^t))" 1 W 2 (t) = (W 2 (t)) t = (WtJI 7 )) -1 V^ 2 (t 7 ) = (W 2 (t 7 )) t (Met 7 )]- 1 ' 

from where each equality follows easily. 

Notice that the first two relations would be equivalent in the Hermitian case. 

Lemma 4.4. For any polyradius r G we have 


□ 


-ia t 


, i r x( z ) • (x( z 1 )) 

I = -FT O 


i r x( z1 )-(x( z )) 


d D z =--— /VA ~ _^ l/v ^~ 77 d u Z} 


(27ti) D Jt d (t) ria =1 z a (27ti) ° Jt d ( r) E[a=1 Z “ 

here d z := d z, d z 2 ... d zq denotes de Lebesgue measure in C D . 

Theorem 4.3. Let c, t g V [U i g (1,2} be two points in the Reinhardt domain of convergence of the Laurent series 
ft(z) of the measure d p and let r c , re G be the corresponding poly-radii. Then, the Baker functions satisfy the 
following bilinear equations for every couple of admissible tunes (t, t 7 ) 

T 1 (z,t)(4 7 J(z- 1 ,t 7 )) t ^ D 

JTD(r c ) 


o 


a 


T D (r c 


O 


T D (r c 


O 


nil z a 


W 2 (z~\t)(^ 2 (2, 

,t 7 )) t 

nLi z « 


M/t (z, t) (^(Z- 1 


nLi z * 


^ 2 (z,t) (^(z" 1 



d z = cp 
d D z = (j) 
d D z = d) 


TNip 1 ] 


T D (tv 


T D (r c 


n°=1 Z a 


d D Z = <±> 


$2( 

z-\t)^* 2 {z, 

t 7 )) 1 


nS=i Zq 


$1 

(z.tJ^JZ- 1 



nS=i Za 


$2( 

Z 1 , t) (#i (z, 

t 7 )) 1 


nLi Zq 




.t')) + 


T D (rL 1 ) 


nS=i Zq 


The MVOLPUT evaluated at different times satisfy 

P(z 


l 0(z,t)(6(z 1 ,t 7 )) t e 

Jt d (r c ) ^ 


d D z = o ®(z 1 ,t)(cb(z, t 7 )) 

[rf 


te 


iD _ 

d z, 
d D z, 
d D z, 


jD _ 

a z. 


,t(z 


ria=l z a JT D (r e ') lla=1 Z “ 


iD _ 

d z, 


o 6(z 1 ,t)(0(z,t 7 )) 

Jt d (r c ) 


, ut e i(z) ft*) 


d D z=A $(z,t)(0{z-'.t')) t e "^ ;i|Z- ' ) d D z, 

tD („— 1 f 


]la=l Z a J T D (r e 


nS=i Za 

(t+t'Hz-’} of--! 


-i + /Ate 


o 0(z,t)(0(z- 1 ,t 7 )) 


(t+t')(z) 


W z ) 


.D _ 

d z = 


TT D (r c ) ' ' EIa=1 z a JT D (r 

I ®(z,t)($(z- 1 ,t 7 )) t d D z = | 

(r c ) ]la=1 Z a JTD (r 4 


= o 0(z~\t)(&(z,t 7 )) 

ttD f 


te 


nS=i Zq 


W z J jD 

a z, 


0(z- 1 ,t)($(z,t 7 )) t ^ ] 

') 1la=l z a 


.D _ 
d Z. 
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Proof. For the first set of equations involving Baker functions just insert the identity according to Lemma 
4.4 matrix in each equality in 4.11 and remember Definition |4.7[ For the second set of equations involving 
MVOLPUT use Proposition 4.8 □ 


4.5. Miwa shifts and vertex operators. For D = 1, see 0, Miwa shifts generate discrete flows from con¬ 
tinuous flows, there you have coherent shifts of the form t —> t 4- [w] or tj —> tj + wVj, j G Z>, where we 
request wGD,zGlD)orwGC\lD),zGC\D; recall that when considered as a perturbation of the measure 
we evaluate z G T. This suggests to consider for each orthant different Miwa shifts. We need to use a 


slightly modified version of Definitions 2.16 and 2.17 and also introduce orthant coherent Miwa shifts 


Definition 4.8. For each subset a G 2* 11 ice introduce 

(1) The corresponding complete integer orthants are given by 


D 


(Z u r : =xz\ 

i=1 


t , Z_, is a, 


Z l := 


Z^_j i G Ccr, 

(2) For any w G (D D ) ff the corresponding orthant Miwa coherent shift [w] ff has as its entries 

„ . , fo, a € IP \ (Z D r, 

«s( z D r- 


Observe that (Z D ) ff = Z D n 
Lemma 4.5. For each cr c Zd and w G (D D ) ff we have 


Y ([w]ff) a z“ = log (1 - Y (-WjZj) 1 - Y_ W i Z 0 
aez D leer ieCw 


uniformly for z G (O d )o-. 

Proof. In the one hand we have 

(4.9) 

and on the other the expansion 


([wlo-ja- 

aez° 


z a = - 


L 


ae 


oc 


— 1 , _ _ V 

log (i - y( w i z i r 1 -X WiZi ) = ~Y ]-aY( wjZjY + X WiZ{ ) 

k=l iger 


i£cr 


iGCo 


ieCo 


□ 


converges uniformly for z G (1L :D ) (J and coincides with the RHS of (4.91. 

Proposition 4.12. For each cr G 2 Zd and each w G (D d ) ct the orthant Miwa coherent shifts in the continuos flows 

tut- [w ]u 

induce, in terms of the Laurent polynomial L w := 1 — Y_ ~~ H w iZi of longitude £(L) = 1 the following 

lea igCff 

Darboux perturbation of the measure 

dp t (0) -> L w (e 10 )dp t (0). 

Proof. From the previous Lemma we have 

exp(t(z) - [w](j(z)) d |x(z) = (1 - ~ X WlZi ) ex pW z )) d h( z ) 

ieo- ieCo- 

and the result follows. □ 
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This Miwa shift can be written as a vertex type operator F^, for each orthant {lP) a and each vector 

wg(D d ) ct 


C : = exp ( - Y- 

ae(Z°) CT 


w“ 9 \ 
|a| 9t a / 


Observe that these vertex operators generate discrete flows or Darboux transformations. Consequently 
they are useful in results like the Miwa expressions of the MVOLPUT, see Theorem |4.1| 


4.6. Asymptotic modules. We will proceed to study families of nonlinear partial differential-difference 
equations involving a single fixed site, say the k-th position, in the lattice and therefore not mixing several 
sites as the Toda type equations do mix, involving near neighbors k — 1, k and k + 1. . 

Since we are interested in partial differential-difference equations in this section we have to consider 
the full deformation matrix, this is, the product of both the continuous and discrete parts that have been 
introduced in previous sections. 

Definition 4.9. For the perturbation parameters t = (t a } agZ D, t a £ C and m £ Z 2D we introduce the following 
complete perturbation 

2D 

W 0 (t,m) := exp ( Y_ t « T «) (j* a ~ qa) m “ • 

aeZ D a=1 


Observe that 


and 


2D 

(Wb(t,m))t = exp( Y_ t-« T a) - qa) m ' 

a£Z D a=1 


9W 0 f = /8W 0 y 
9t a V 9t a ) 


T Q WJ = (T a W 0 ) t . 


Definition 4.10. Given two semi-infinite matrices R] ft, m) and R 2 (t, m) we say that 

• R] (t) £ IWo z/Ri (t) (Wo(t, m)) 1 is a block strictly loiver triangular matrix. 

• Ri(t) £ (IWo)t z/Ri (t) ((Wo(t, m))^) 1 is a block strictly lower triangular matrix. 

• R 2 (t) £ u if it is a block upper triangular matrix. 

The sets [Wo and u are known as asymptotic modules. 


Then, we can state the following congruences l66l or asymptotic module Il54l style results 


Proposition 4.13 (Asymptotic modules). (1) Given two semi-infinite matrices Ri(t, m) and R 2 (t, m) such 
that 

• R^tj-m) £ [W 0 (t,m), 

• R 2 (t,m) £ u, 

• Ri (t,m)G = R2(t,m). 

Then, Rj (t, m) = 0 and R 2 (t, m) = 0. 

(2) Given two semi-infinite matrices Ri (t, m) and R 2 (t, m) such that 

• Ri(t,m) £ (lWo(t,m)) t , 

• R 2 (t,m) £ u, 

• Ri(t,m)G^ = R 2 (t, m). 

Then, Ri (t, m) = 0 and R 2 (t, m) = 0. 


Proof. Observe that 

R 2 (t, m) = Ri (t,m)G = Ri (t, m)(Wi (t,m)) 1 Wi (t, m)G = Ri (t,m) (Wi (t,m)) '^(t, m), 
R 2 (t, m) ^H(t, m) ((§(t, m)) _1 ) f ) = R, (t, m) (W 0 (t, m)) _1 (S(t, m)) _1 , 
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and also that 

R 2 (t, m) = Ri (t,m)G^ = Ri (t, m) (Wi (t,m)) (t,m)G^ = Ri (t,m) (YVi (t,m)) 1 W 2 (t, m), 

Ri(t, m) (hl(t, m) _1 S(t, m)^ f = R, (t, m) (W 0 (t, m) t )~ 1 (§(t, m)) _1 , 

and, as in every RHS we have a strictly lower triangular matrix while in every LHS we have an upper 
triangular matrix, both sides must vanish and the result follows. □ 

We use the congruence notation 

Definition 4.11. When A — B G IWq (g (IWo)t) we write A = B + IWo (A = B + (IWo)t) and if A — B G u we 
write A = B + u. 

4.7. KP flows. We now study, using the asymptotic module technique, the first and second order KP flows 
and their compatibility conditions resulting in nonlinear partial differential-difference equations for the 
coefficients of the MVOLPUT. These are onsite equations not mixing different values of k. Finally, we 
consider the third order flows. 


Definition 4.12. 


where 


The longitude one perturbation parameters are 

f[i] :=(t-b--->t_D,t D ,...,t 1 ), 


t Q t sgn ( Q ) e|a| , 

We also consider the following derivatives 


a G {1,..., D}. 


3 n := 


9n := 


3n a ■— y Ll Qi b3b, 


3n a •— y Ll a b3b, 


where n a G C 2D , a G {±1,..., ±D}. 
Proposition 4.14. The following relations hold true 
3n a W! =[(I + (3)(7t Q )>]Wo + [W 0 , 
3n Q W T =[(I + 0)(ftJ>]W ( 5 + ([W o )t, 
Coinciding both equations in the Hermitian case. 
Proof. Firstly, we realize that 

3n a W 0 =7f Q W 0 , 


TaW! =[(I + T Q (3)(Ti a )>-q Q ]Wo + [W 0 , 
T q V^ —[(I + T a 0)(ft a )> - q Q ]WJ + (IWo)t. 




Secondly, we remember the definitions of Wi and Wi and work there using the congruence techniques 

3n a W T =[3 na S + STf Q ]W 0 

= [(I+(3)(|f a )>]Wo + lWo, 

T q W! =(TaS)(Tf Q -q a )W 0 

-[(I + T a (3)(Tt a )>-q Q ]Wo + lW 0 , 

3n a V^i =[3 flQ $ + ^ a ]Wj 

=[(i+^)(ft a )>]w 0 t + aw 0 )t, 

T Q ^ = (T a §)(fl a -tJa)Wj 

= [(! + - q Q ]wJ + (lW 0 )t. 

□ 
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Proposition 4.15. Both Baker functions {^ 1 ,^ 2 } are solutions of the following partial difference-differential linear 
systems 

(3 na -T a )¥ = [q Q -A a (3(it Q )>]V. 

While both Baker functions {^1 ,^ 2 } do solve 

(3n a - T a )W = [q Q - 

Only one set of these equations is relevant in the Hermitian case. 

Proof The results in Prop osition |4.14 1 can be rewritten as it follows 

(3n Q -T Q -[q a -A Q |3(TfJ>])W 1 e IW 0 , (3 fta - T a - [q a - A € (lWo) t . 

Then, it is straightforward to see that 

(3n Q - T Q - [q Q - A Q (3()f a )>])V^2 € u, (3 fla - T a - [q„ - A a 0(ftj>])W 2 G u. 

Finally, just recall the definitions of the Baker functions and apply Proposition |4d3] □ 

The compatibility of the linear systems satisfied by the Baker functions imply 
Theorem 4.4. The following nonlinear discrete equations for the first subdiagonal matrices |3 and (3 hold 

(4.10) A b [3n a (3 + (A a (3)(q a + (rtJ>(3)] {yi h )> = A a [3 nb |3 + (A b (3)(q b + (* b )>p)] (tf Q )>, 

(4.11) A b [3 ft J + (Aj)(q a + (ftj>0)] (ft b )> = A a [d n J + (A b $)(q b + (ft b )>0)] (ft a )>. 

Definition 4.13. (1) Longitude 2 times are couples (a, b) c (—1,..., D, D,..., 1} where we assume a + b f 0 


t(a,b) : = tsgn(a)e| a |+sgn(b)e| bj ) 

with corresponding derivatives denoted by 

3 
3 r 


J (Q,b) 


3ti 


a,b) 


(2) We introduce the diagonal matrices 

(4.12) V Qjb :=3 a (3 (T b ) > , (V a>b( )=3 a (3 [lc ] (Yb)[ 1c _i] i [ lc ] > h a , b := — V a , b — V b , a , 

(4.13) ^a,b : = 3-a0 (T b ) > , (^a,b)[k] =3-a0[k] Pf b )[k-1],[k] > ^a,b := ~ ^a,b ~ ^b,a- 

We have excluded the case a + b = 0 because in that situation what we obtain is t( a _ a j = to. Observe 
that 


3¥ t 
3 to 


= V* 


^i = , p . 

3t 0 


with 1 = 1,2. 

Proposition 4.16. Both Baker functions d ; i and ^2 ar e solutions of 

dV d 2 V 


(4.14) 


3t(a,b) 3t Q 3t b 


v = 


d 2 v 


3t a 3t_ a 


U a , b T, a + b/0, 


+ U a ,_ a ¥, 


while Baker functions and V 2 solve 


3^ 


d 2 W 


(4.15) 

Proof. See Appendix |C.11| 


3t ( _ a _b) 3t_ a 3t_ b 
d 2 w 


+ ti a)b ¥, a + b/0, 


V = 


3t_ a 3t a 


+ &a,-aV. 


□ 
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Observe that for a = b ( |4.14| ) reads 

3%] 3 2 %] 


9t 


( 2 ) 


dt 2 a 


and (4.151 


9¥ 


[k] 


d 2 Wi 


[k] 


3tS 


9t 2 , 


+ (U a ) H %], 


+ 0X)[k]%], 


+ ra . , 

L a • L (a,ap 


I(2) 

L —a • L (—a,—a)) 


Ua=- 2Vn 


da=-2Q a , a 


which are time dependent one-dimensional Schrodinger type equations for the square matrices Ty, the 
wave functions, and potential the square matrices (U a )[ k ], (H a )[ k ]. Moreover, multidimensional matrix 
Schrodinger equations appear if we look to other directions, thus given (ai,..., ay) C {±1,..., ±D}, we 
can look at the second order time flows generated by := 


at 


(“1 ,a,) 


+ ••• + 


at 


( Q d> a d)' 


A 

at 


at 


(-“i ,-<m) 


+ ••• + 


yy— -- to get in terms of the following d-dimensional two possible nabla operators V := (gT—,..., g^— 




V := 


at_a, ’ ’ ’ ■ ’ 3t_ 


Laplacians A := V 2 = 


+ ••• + 


-, A := V 2 = 


9t a, 1 1 v 91^, ' ' Mu Qd 


+ ■•■ + 


and matrix 


[k] 


= &y K + n rkl ^ 


[k] ” [k] • 


potentials U := U Q1)Q1 H-h U ad)Qd = 2V((3) • T, 

9Vr kl 9^ 

-gf = A %] + U M%] -g r 

Where the equation on the LHS suits M ; i and ^2 while the one in the RHS suits $1 and MG. 

We see that again only the k-th site of the lattice is involved in these linear equations and, consequently, 
its compatibility will lead to equations for the coefficients evaluated at that site. These nonlinear equations 
for which (3 [y or |3 [y are a solution are 


Theorem 4.5. The following nonlinear partial differential equations 

9 ( c,d) (3a(3(T b )> + 9 b |3(T a )>) - 9 (a>b) (9 c (3(T d )> + 9 d (3(T c )>) 

= 9 a 9 b (9c|3(T d )> + 9 d (3(T c )>) - 9 c 9 d (9a(3(T b ) > + 9 b (3(T a )>) 

(4.16) + (9 b 9 c |3)((T d ) > (3(T a )> - (T a ) > |3(T d ) > ) + (9 b 9 d |3) ((T c ) > (3(T a )> - (T a ) > (3(T c ) > ) 

+ (3 Q 3 C |3)((T d )>(3(T b )> - (T b )>(3(T d )>) + (9 a 9 d |3)((T c )>(3(T b ) > - (T b ) > P(T c )>) 
+ [3 a (3(T b )> + 9 b (3(T a )>, 9 c (3(T d )> + 9 d (3(T c )>] 


3(- c ,-d) (3-a0(T b ) > + 9_ b j3(r Q )>) - 9 ( ^ b) (9_J(T d )> + 9_ d 0(T c )>) 

= 9- a 9_ b (9_ c 0(T d )> + 9- d $(T c )>) - 9_ c 9_ d (9_ a $(T b )> + 9^ b 0(T Q )>) 

(4.17) + (9_ b 9_ c (3) ((Y d ) > 0(T a )> - (T a )>^(T d )>) + (3_ b 8_d0)((r c )J(r Q )> - (T a ) > ^(T c ) > ) 

+ (9_ Q 9- C $)((Y, d ) > 0(T b )> - (T b )>0(T d ) > ) + (9_a9- d $)((Y c )>$(Y b )> - (T b ) > 0(T c ) > ) 
+ [3-J(r b )> + 9- b 0(T a )>, 9- c 0(T d )> + 9_ d 0(T c )>] 
are satisfied for a, b, c, d e (±1,..., ±D}. 


Observe that this equations decouple giving for each k the same equation ( |4.16| ,( [4TT| ) up to the replace¬ 
ments (3 -> (3 [k] ,0 -> j3 [k] and T A -> (T A ) [k _i ])[k] , k = 1,2,... and A = a,b,c, d. 

Proof. Let us start with the first equation. If we denote 

(4.18) L Qjb := 9 Q 9 b + U ajb , a, b e (±1,..., ±D}, 

( 4.14} reads 9 a , b (Wi) = L Qjb (Wi). The compatibility conditions for this linear system are 

^3(a,b)(h c , d ) — 3( c , d )(L a , b ) + [L C)d ,L Q)b ]^(Wi) = 0 


and consequently 

(4.19) 


Ra, b ,c,d(Wi) — 0, 


a,b,c, d e {±1,...,±D}, 
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where 


Ra,b,c,d : = (3bU Cj d)3 a + (3 Q U C) d)3b — (3dU Q) b)3 c — (3 c U Q) b)3d — 3( a ,b) (hl c? d) + 3(c,d)(U a , b ) 

+ [U c ,d)U Q) b] — 3 c 3dU a> b + 9 a 3 b U C) d. 

Therefore, Wj satisfies an equation of the form 

(^B j 3 j +A)(W 1 ) = 0. 
jei 

for a given subset of indexes I C {±1,..., ±D}. Now, recalling that 

3aW, = ((T q )> + |3(T a )>)Wo + [Wo 

leads to 

0= (^B j 3 j +A)(W 1 ) = (^Bj(T j ) > + ^DB j (3(T j ) > + A)Wo + [Wo 
jei jei jei 

Observe that Y_ Bj (3(Tj)> + A is a diagonal matrix. The sum Y_ Bj (Tj)> is a linear combination of terms like 

jei jei 

3a(3fT Q )>(T b )>, which either are zero everywhere but possibly on the first superdiagonal a + b / 0 or are 
identically zero for a + b = 0, since fT a )>(T a ) > = 0. Therefore, each of them can be set equal to zero and 
from ( |4.19| ) we find, in the first place, that 

(3 b U c , d )(T Q )> + (3 Q U c ,d)(T b )> - (3 d U Q) b)(T c )> - (3 c U Q)b )(Td)> = 0, 

which is identically satisfied because of ( 4.12) . In the second place, we get the following nonlinear equation 

(3bU c ,d)(3(T Q )> + (3 Q U c ,d)|3(T b )> - (3 d U Q , b )|3(y c )> - (3 c U Q)b )|3(r d )> + 3 (Q)b) U C)d - 3 (C)d) U Q , b 

+ [Uc,d> hl Q)b ] — 3 c 3 d U a)b + 3 a 3bU Cid = 0, 

and recalling ( |4.12) we get the desired result. 

The proof for the second equation is exactly the same with the repacements 3 Q —> 0_ Q , (3 —> U Qib —> Q Q)b , 
[Wo —> (IWo) t . ’ b 

Definition 4.14. (1) Longitude 3 tunes are triples (a, b, c) C {±1,..., ±D} where zve assume a + b / 0, a + 

c/ 0 , b + c 7 ^ 0 

f(a,b,c) := tsgn(a)e| a |+sgn(b)e ib |+sgn(c)e| c |> 

with corresponding derivatives denoted by 


(a,b,c) 


3t, 


a,b,c) 


(2) We introduce the diagonal matrices 

^a,b,c •“ diag((V aj b,c)[ 0 ]) (^a,b,c)[l]) (^a,b,c)[2]> • • • ) 

with 


V a ,b,c :=3a(3[(3, (T b ) > ](T c ) > , 

(V Q)b ,c)[k] =3aP[k](p[k-1](^'b) [1c _2],[k-1] ~ (^ b )[k- 1 ],[k]Pw) ( T c) [k ^ 1])[k] 

We have excluded the case a + b = 0, a + c / 0, b + c / 0 because in those cases what we get is 
t( a ,b,c) = tc, W, t Q , respectively. 

Proposition 4.17. Observe that 

Va,b,c =3a|3 (2) (T b )>(T c )> - 3 a |3(T b ) > (3(T c ) > , 

(V a ,b,c)[k] =3a(3[ k] ) (^b) [k _ 2 ],[Tc_l] _ 3 a (3[ k ](Tb) [k _ 1 ] )[k ]|3[k](T c )[ lc _ 1 ] j [ lc] . 
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Notice that there is a similar definition ^ a> b,c given by the modifications |3 —> 0 and 3 Q —> 3 Q ; for the 
sake of simplicity, we will just mention it after every result instead. Observe also that to obtain the second 


equalities we have used ((3 a S)S , = 0, which follows from Lemma 


4.3l that leads to 


® Q P[kl — ^ Q P[kl P[k—1]) 


o(2) ._ O 

P[k] — b [k][k-2]* 


( 2 ) 

We remark that (V a> b )C )[ k ] depends on |3[ k ] and (3^ only, coefficients of the MVOLPUT for the second 

and third higher lenght monomials, O m (z) = X[k] (z) + P[k]X[k-l] ( z ) + p|yX[k- 2 ] (z) H-h (3^. If we insist 

in using only the second higher longitude coefficient and not the third one there is a price that must be 
paid, two consecutive Laurent polynomials 0[ k ] and ®[ k -i] will be involved —as we require both of |3 k j 
and |3[ k ^i]. Then 


Proposition 4.18. The Baker functions Ti and $2 are both solutions of the third order linear differential equations 
3(a,b,c)^ =3a3 b 3 c ¥ - (1 - 6 0 ,b+c)V Q , b 3c^ - (1 - 6o,a+c)V b)C 3 a ^ - (1 - 5o,a + b)V c , a 3 b ¥ 

— (9aV b) c + 3 c V a(b + 3 b V Cja + (1 — 6o,b+c)V a)b ,c + (1 — ^0,a+c)Xb,c,a + (1 — So ) a+b)V C) a )b )' V l / . 
and a similar equation holds for ^ and ^2 a ft er the replacements V —> V, (3 —> $ and 3 a —> 3^. a . 


Proof. See Appendix |C.12| □ 

Notice that as we anticipated in the study of the second order flows there is a special feature to be 
considered with rescpect to what we would expect from the real case of MVOPR Ifl3l . As one can see, the 
equation takes different shapes depending on the values of a, b, c. 

• If |e Q + e b + e c | = 3 then none of the 5's is present and the result coincides with the one we had for 
MVOPR. 

• If |e Q + e b + e c | = 1 then observe that in the LHS of the equation we would be deriving with 
respect to a first order flow (instead of a third order flow). Also at least one, and at most two of the 
following three equations must hold b + c = 0, a + c = 0 and a + b = 0. As a result at least one 
and at most two of the 5's give a nonzero contribution and two of four of the terms in the equation 
cancel. 


Appendices 


A. OLPUC 


In 0, following a CMV approach, orthogonal Laurent polynomials in the unit circle (OLPUC) where 
discussed in detail from the same viewpoint as in this paper. Let us investigate how OLPUC are recov¬ 
ered as MVOLPUT when D = 1 (MVOLPUTd = i). For that aim we recall now some basic facts about the 
construction presented in 0. The moment matrix considered in 0 can be expressed as 


G := 


X d PX f > 

T 


( 1 \ 

^-1 



z^ 1 


( G o,o 

Go,i 

G 0 ,(k.-1) \ 


z 

G m := 

Gyo 

Gi,i 

Gi,(k-i) 

X(z):= 

z -2 



r 2 






z 


\G(k-i),o 

G (k-1),l 

•• G (lc _ 1)>(k _i)/ 


Where x is the vector of monomials ordered according to the CMV proposal. Notice that as we assume here 
that the measure is real, d q = d jl, we have that the moment matrix is Hermitian, G = G ' ; moreover, we 
also assume that we are dealing with a positive definite Borel measure, i.e., every minor of the associated 
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moment matrix is positive, this is det > 0 for all k <G {0,1,2,...}. For such cases the Gauss-Borel 
factorization of the moment matrix leads to 


(A.l) 


G := s 1 h(s 

/I 0 0 


s := 


«i 

* 


1 0 

CX2 1 

* 


«2k+1 

* 


1 0 

«2k+2 1 


V 


h := diag(h 0 , hi,..., h 2k+ i h 2k+2 , • • • ’ 




with the OLPUC given by 

<p(z) := s X (z), cp (2k+1) = z~ (k+1] + • • • + n 2k+1 z k , 

and corresponding Szego polynomials expressed as follows 

'd c +L n (2k+1) _ p* 


<P 


(2k+2) _ 


— « 2 k+ 2 Z 


•—(k+ 1 ) 


+ • • • + z 1 


k+1 


<P 


(z) = P 2 k+i(z), 


zk +1 Cp ( 2k + 2)(z)=P 2k+2( z), 


with reciprocal polynomials given by P£(z) = z k P|Jz 1 ) and Verblunsky coefficients defined as a k := P k (0). 
In the theory we have also p k := 

However, the MVOLPUTd=i is not exactly the one just described for OLPUC 0. Despite we are dealing 
with the very same moment matrix G, now the Gauss-Borel factorization is a 2 x 2 block (but for the first 
left upper corner) factorization: 


(A.2) 

with 


/I 0 0 

Pm h o 

* p[ 2 ] 1 2 


S := 


G = S~ 1 H(S“ 1 ) t , 
\ 


P[k+1] h 


V 


/H [0 ] 0 

0 H m 
0 


H := 




0 
0 

0 h [2] 


H 


[k+1] 


V 




being all the entries of these matrices, except from S [o] [o] = 1, H- 0 , which are real numbers, 2x2 complex 
matrices. The Laurent orthogonal polynomials are encoded in the semi-infinite vector 


4>(z) := Sx(z), 


where we can write 


x(z):= 


/ 1 \ 

X[l](z) 

X[ 2 ] (z) 


X[k+l](z) 


V 


X[k+i] : = 


,-(k+ 1 ) 

^k +1 


/ 
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and 


k+1 

4)[k+i]( z ) : = S[ k +ii[nxm( z ) = 

1=0 

To properly connect both scenarios we need a further object: 


f cj)* 2 k+ 1 Hz)\ 
\^(J)t 2k + 2 )(z) J 


/ z -(k+ 1 ) _|- + * z k\ 

V *z k -I -b z k+1 J ' 


Definition A.l. The derivative reciprocal Verblusnky coefficients are defined by 

d , i (0) 


^ 2 k+i := 


2 k+1 1 

dz 


Then, we have the following 

Proposition A.l. When the Gauss-Borel factorization ( |A.1[ ) can be performed then ( A.2| also exists with the block 
entries ofS , H expressed in terms of the entries of s, h as follozvs 

1 «2k+2^ Q _O^k+l^ 

«2k+2 1 J ’ 


H[k+1] = —h2k+i 


P[k+1] = 


\ a 2k+1 ^2k+1 J ’ 

where in the last equation we take k ^ 1 and for k = 1 we have (3[i] = ^ “J j. Moreover, 

Ai = ai A2k+3 — A2k+1 = «2k+1«2k+2 + «-2k+2«2k+3) 

it also holds that 


Pk — 1 ^k^k- 

The MVOLPUTd=-\, {(^(z)}^, the OLPUC, {(pkjz)}^, an d ^e Szego polynomials mid reciprocal polynomials, 
{P k (z), Pk(z)}^To' are re l a t e d as follows 


(j)f 2k+ b( z ) = (pf 2k+1 *(z) 


ff 2k+2 \z) = — CX2k+2<P (2k+1) ( z ) + Cp ,2k+2) (z) 


A +1 o 


=(\ 


H-' 




Proof. To prove the results we must start removing from s its block diagonal. We introduce 

/I 0 0 0 0 \ 

0 10 0 0 

0 1 x 2 1 0 0 

0 0 0 10 

0 0 0 04 1 

n := 


1 0 

«2k+2 1 


V 




If we do so in the factorization we end up with 

G = s~ 1 nn~ 1 h(n~ 1 ) t n t (s~’ 1 ) t = (n^s) -1 (n^h^ -1 ) t )[(n" 1 s)- 1 ] t , 

so that 

S=u~ 1 s, H = n^'hjn -1 )t. 

From the first expression the relations between cf> and cp follow and from the second one we get 

H = ( ^ 2k+1 h 2 k+1(X2k+2 \ 

[k+]] \«-2k+2h2k+1 h 2 k+2 + «2k+2h.2k+l «2k+2/ ’ 
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But due to the persymmetry of H[ k+ i] the elements in the main diagonal must coincide and therefore 

h2k+1 = h-2k+2 + «2k+2h2k+1 &2k+2> 

from where we deduce P 2 k +2 := nyfy = 1 ~ « 2 k+ 2 « 2 k+ 2 - In order to get the expressions relating the 4 ) to the 
Szego polynomials we only have to take the relations between cj) and (p and remember both the expressions 
relating the <p to the Szego polynomials and also the recursion relation for the last ones 

Pk(z) = zPk -1 (z) + CXkPk -1 ( z )- 


These relations lead to expressions for the entries of Ppk+ip From the two alternatives of writing the Jacobi 
operator we get the relation between a and A 


OC 2k+1 


J[k+1],[k+1] a 2 k+3 A 2k +1 — A 2k +3 

_ h2k+l / A 2 ic + i— A 2 ic + 3 — a 2 ic + l a 2 k+2—«2k+2a 2 k+3 ~«2k+2(A 2 k+1 ~ A 2 k+3—«2k+2«2k+3)+&2k+1 

h"’k+'’ \ [ a 2k+2(A 2 k+l— A2k+3)—«2k+3]~ 0 C 2 k+ 2 a 2 k+ia 2 k +2 [« 2 k+2(A 2 k+1 — A 2 k+3)“«2k+3]&2k+2+<*2k+2a2k+1 J 


While from the recursion relation for (J), which can be expressed as follows 


z4>[k+i] = 


h2k+l 


1 


-a 2 k 


h2k V a 2k+2 OC2k+2«-2k 


4* [k] + 


0 


a2k+1 


— «2k+3 <*2k+1 «2k+2 + a2k+2«2k+3 


4>[k+1] 


0 0 

0 1 


4* [k+2] > 


we get 


z4*2k+1 (z) = " k+1 [4>2k —1 _ a 2k4 ) 2k] + a2k+l4 ) 2k+2- 

h 2 k 


If we compare the coefficients in z k we get 


, h 2 k+i _ 

1 = —r -P <*2k+1 a 2k+1 

n- 2 k 


which, recalling that p 2 k+i : = h ^ k+1 gives the expressions for the odd p. 


h-2k 


□ 


Notice that for the previous result to hold we have assumed that the Gauss-Borel factorization ( A.1J 
can be performed, which implies the block Gauss-Borel factorization ( A.2| l. The OLPUC factorization 
( A.lJ of 0 was ensured whenever all the standard principal minors do not vanish, det G k / 0 for all 
k E {0,1,2,...}; however, for the MVOLPUTd=i construction we request a less demanding condition, 
now we require only the block principal minors not to vanish, hence it only requests det G 2k| F ^ 0 for 
all k E (0,1,2,...}, and nothing to the even minors, which in fact could vanish. Whenever all the odd 
principal minors are non zero and some even principal minors do vanish we have MVOLPUTd=i but not 
OLPUC. Let us explore this situation 


Proposition A.2. Let us consider a measure such that det G t2k+ F / 0, k E {1 ,2,...}, det G^ 2k+2 ^ / 0 Vk / l and 
det G [2l+2] = 0. Then 


and 


H [k+1] - 


/ h 2 k+i 
\ a 2k+2h.2k+1 


h-2k+l «2k+2\ 
h-2k+1 J 


H [l+1] - 


/H 2 i+1,21+1 

\Fl2l+2,2l+1 


H21+1,21+2 
^21+2,21+2 


(_0 H21+A 

V H 2l+1 0 J 


Vk^l, 


H 2 i+i := (ct> (21+1) ,4> (2v+2 )). 
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Proof. From the imposed conditions on the moment matrix we see that Vk < 21 + 2 we have that det g ^ >0 
and therefore det H[ k ] > 0. But det g 21 2 = 0 and therefore 


H [0 

0 

0 


0 

H[i] 


0 

0 


0 H 


[ 2 ] 


Hr, 


Hh+1,21+1 


L 

det H m j H 21 .+ 1 , 21+1 = 0. 


k=0 


Since f7k_ 0 det 1 l[ k ] = det g 2L+1 ' > 0 the only possibility is that H 21 + 1 , 21+1 = 0- Then using the persym- 
metry of the blocks of H it is straightforward to get that 1 121 , 2 , 21+2 = H 21 + 1 , 21+1 = 0 an d that 1 121 , 1 , 21+2 = 
Hu 1 - 2 , 21-1 := H 21 + 1 . Notice that since again det G 2l+3 l > 0 we have that det 1 lf lfl f 0 and therefore 
H 21 + 1 V 0 . □ 


We have 


_ f h2k+i h2k+ia2k+2\ _ ( 1 0\ /h2k+i 0 \ ( 1 oV 

lk+1] V a 2k+2h2k+l h-2k+1 ) \ a 2k+2 1/ \ 0 h 2 k+l ) V a 2k+2 1 ) 

B. Persymmetric and spectral matrices for D = 1,2 
For D = 1 the per symmetry matrix is For D = 1 we have 


( 

1 

0 

0 

0 

0 

0 

0 

0 

0 

... \ 


0 

0 

1 

0 

0 

0 

0 

0 

0 



0 

1 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

1 

0 

0 

0 

0 



0 

0 

0 

1 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

1 

0 

0 



0 

0 

0 

0 

0 

1 

0 

0 

1 



0 

0 

0 

0 

0 

0 

0 

0 

1 



0 

0 

0 

0 

0 

0 

0 

1 

0 


V 





/ 


and D = 2 the partial persymmetries matrices are 


/ 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

... \ 


0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 



0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 



0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 


V: 



) 















1 000000000000 
II 0 0 0 0 0 0 0 0 0 0 0 
0001000000000 
0010000000000 
0000100000000 
0000010000000 
0000000100000 
0000001000000 
0000000001000 
0000000010000 
0000000000010 
0000000000100 
0000000000001 


For D = 1 we have the spectral matrix 


0 

0 

1 

0 

0 

0 

0 

0 

0 


1 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

1 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

1 

0 

0 


0 

0 

0 

1 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

1 


0 

0 

0 

0 

0 

1 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 



As an example, we show the first blocks for the D = 2 case 


0000100000000 

Xoooooooooooo 

0000000000100 

0000000000010 

0000000000001 

“(TlOOOOOOOOOOO 

0010000000000 

0001000000000 

0000000000000 

0000000000000 

0000000000000 

0000000000000 

0000000000000 
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( 0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 


1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


V : 





C. Proofs 

C.l. Proof of Proposition |2^4j 

Proof. (1) Assuming det A f 0 for any block matrix M = ( + f] ) we can write in terms of Schur com¬ 
plements 


M = 


I 0\ (A 0 
l -1 


CA-' I) VO M/A J VO I 


I A -1 B 


Thus, as det G lkl / 0 Vk = 0,1,..., we can write 
G [l+1] = 



0 ^ 

( G 111 

0 ^ 

( I 111 

«)[l],[l—1])T 

v [l],[l—1] 

%]l J 

l 0 

G 11+11 /G 111 J 

l 0 

%]l 


where 


v [ll,[l 11 := (v [l])[0] v ra>[1 ] ... 

Applying the same factorization to G' L we get 



/ ][[l-1] 

0 0 \ 

( G 11 - 11 

0 0 

G 11+11 = 

r [1-1] [1-2] 

%-1]l 0 I 

0 

G 111 /G 11-11 0 


^ s [l] [1—2] 

■t[l][l-1] Il[l]l J 

l 0 

0 G 11+11 /G 111 



/ ][[l-1] 

^[l-1][l-2]jT ^[l][l-2]^T 

X 

0 

I i[i-i]i fti],[i-i]) T 


l 0 

0 %]l 


Here the zeros indicate zero rectangular matrices of different sizes. Finally, the iteration of these 
factorizations leads to 


Aral 0 ••• 0 \ 


G 11+11 = 


* In 


V 


A ion o • • • o \ 


T 


diag(G [11 /G [0] , G 12 ] /G m ,... ,G LL+I] /G 


[i+i] /(-'[ i]i 


%]|/ 


* In 


V 


0 

* %][/ 


Since this would have been valid for any l it would also hold for the direct limit lim G 111 . 
Consequently, = G ll+11 /G 111 is a Schur complement and det G ll+11 = nlc=o det H- k , f 0. 



































Thus, 


G f = (o x(z(0))dn(0)(x(z(0))) 


T D 


= o x(z(0))dp.(-0)(x(z(9))) t - 
Jt d 


G 


G f 



X(z(0))[d^0) - d p-(—0)] (x(z(0))) f . 


Consequently, for real measures the moment matrix is Hermitian. Conversely, if G is Hermitian 
then d |i(0) — d pi—0 ) has all its moments equal to zero, and therefore can be identified with zero. 

Now, as the factorization is unique from § _1 Ht (S -1 )^ = G^ = G = S -1 H(§ -1 )^ we deduce that 
§ = S and Ht = H. 

(3) The truncated moment matrix G k ' satisfies for any truncated vector v |lc = (v[o], Vpj,..., ) T <E 

Iffll 


k-1 


(v [k] ) t G [k] vM = _^() (v [i] ) t x [ i](z(0))dp(0(x[ j ](z(0))) t v [ j ] 


1 - 0=0 


= o 


T D 


k-1 


T D 


Y ( v [i]) f X[i](z(0)) 


i=0 


d M-(0), 


and as the measure is definite positive we conclude that the truncation is a definite positive matrix. 
Conversely, if every truncation is definite positive then |fic(0)| 2 d pt(0) > 0, for truncated Fourier 
series f^. Consequently, the measure is definite positive. 


□ 


C.2. Proof of Proposition 2.10 


Proof. From 


C=(G§t H - 1 )t x (z) 

=(H- 1 )t§Gt x (z) 


we conclude 


m t c)«= y $«,«, x G a 2 , ai z “ 2 

0<|«i|<|a| a 2 SZ D 


=(2tt) d Y_ §<*,«, Y c «i-«2 Z “ 2 

O^loti |<|a| a 2 SZ D 

=(2tc) d Y_ $<*,«, z ai Y e N z ^ 

0^|ai |<|a| a.' 2 &lP 

= (27T) D $ a (z)jI(Z~ 1 ). 


For the second relation we deduce 


<?=H _1 SG X (z), 
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and therefore conclude 


£_ S a , Kl G ai , a2 z«2 

0<|«i K|a| <x 2 eZ D 

= (2tt) D Y _ S «,«i Y c -«i+«2 Z “ 2 

0<|«i ZM a 2 eZ° 

= (27t) D Y_ S a.a^'* 1 X. C «i Z “ 2 
0<|«i |<|a| ct'eZ D 

=(27r) D 4) a (z)^(z). 


□ 


C.3. Proof of Proposition |2.16[ 

Proof. We need to show that for each multi-index a G Z D there exists a unique cr C Zq such that a G (Z D ) ff . 
To prove it we look at the components of a and construct three associated subsets, & ' C Zd for the strictly 
positive/negative entries and a 0 C Zq that collects the components with zero entries. Consequently, we 
have that C a and cr > C Co, and Z^ = U u 0 U o > with cr < n a 0 = a < n a > = a 0 n o > = 0. To 
construct our set a we need to complete ct < with elements in a 0 , but which? We split cr° = a 0, U cr 0,+ , with 
cr 0,=F C a 0 and cr 0, ‘ n cr 0,+ = 0 in the following manner. Take an element in i G cr ' and check whether or not 
i — 1, i — 2, — p G cr° and i — p — 1 0 a 0 , if it happens add the elements i — 1,..., i — pto cr°> - , after 

completing this procedure Vi G cr < we have constructed a 0, and we are ready to define cr 0, := a 0 \ a 0, . 

D 

Then, we declare ct := ct < U cr 0 ’ - so that Ccr = cr’ LJ cr 0,+ . By construction (Z D ) CT = X Z\ where each Z\ is as 

i=l 

follows, if i G cr 0, then i G (a \ 3a) so that Z l = Z _ = Z> U {0}. The set cr 0,+ is built up of strings of zeroes 
lying at the left of an element in a > , hence cr 0,+ C (Ccr \ 3Ccr) and for i G cr 0,+ we have Z b = Z + = Z> U {0}. 
For the elements i G cr < we could have Z t = Z. for i G cr \ 3cr or Z L = Z<, otherwise; a similar discussion 
holds when i G cr > . □ 


C.4. Proof of Proposition |2.20| 

Proof. First, we prove that the set of spectral matrices {T Q }^ =1 is an Abelian set. In the one hand, a straight¬ 
forward calculation leads to the relations 


I Ml 


((TaT b ) [k+1])[k ^]) y - ^( T a) w (k+l) te( C , c)(Tb) (x (k, )W ( , c- 


r=1 

= 8 a' k+1) +e Q+ e b ,a ) ! k - 1 >’ 
I Ml 


(CY , a Tb)[ k+ 1] ) [ k+ 3])y = L (Ya) a (ic+i) a | k + 2 ) (Tb) a f+ 2 'i !a f k + 3 ) 


T=1 

= 6 a[ k+1) +e Q+ e b ,a' k+3 >- 


On the other hand 


I Ml IPc+2]| 

Xf T ^af k+ 1 ) ,a' k)(Tb ^' k) ! a[ k+1) + Y a f\af +2) (Tb) a[ k+ 2 ) ,a ' k+1 

r=l r=l ' 

IM 

= Z 6 „< 


r=1 


x| k+1 ) +e a ,a| k) 8 a| k) +e b ,a! k+ ' 1 


|[k+2]| 

+ Y 5 a! k+1 > 


r=1 


„( k +2)5 (k+2) (k+ll- 

+ 6a,0C r £X r 


In the RHS the first term involves 5 ( k +o ( k +n which is not zero only when the a-th component of 

the multi-index aj k+1) is negative and the b-th component of the multi-index cxj k ' 1 is positive; not only 
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but in the second term of the RHS for non zero terms we need the a-th of a[ k+11 to be positive and the b-th 
component of a! k+1 * to be negative. Hence, we are left with 


I Ml 

((TaTb^+iyk+u^ = ^(T a ) a ( k+1))a ( k )( T b ) a ( k))a ( k+ i) 
r=1 


— S (k+l) „Cle+1) 

a, +e Q +eb ) tt i 


Second, we show that Y a is an orthogonal matrix; i.e., Y a Y a 
the Y's we deduce that 


[k+2]| 

Y ( T ^af\af +2) (Tb) a' k+2) ,«' k+,) 

r=1 


Y a Y a = I. From the very Definition 


2.20 


of 


0^1 )„W a (k+ 1 ) 

w i 

rTi 


( T aJ a ('<+l) )a ( k) 


- 6 a' k) -e a ,«; k+1 )> 
= 6 a{ k+ ’D ea)a W- 


We compute Y a Y a as we did with Y a Y h and we get the identity. Finally, we observe that 

(haY^a^k) (k+O = \ ia [ k G„, , „(k+1) = ^(k)^ T ^ „(k+l) 


x{ +e b ,I a «j 


+I a e b ,«i 


But as I a e b = e b whenever a / b and I Q e a = —e Q we get the desired result. 


□ 


C.5. Proof of Proposition 2.23 

Proof. (1) From the string equation ( |2.13| ) and the Gauss-Borel factorization 

Y a S^ 1 H(§- 1 ) t = S -1 H(§ -1 )^Y a 

and hence 

SY a S- 1 H = H(§- 1 ) t Y a (§) t , 

and the result follows. 


we get 


(2) From the Definition 2.22 we deduce that J a and J a have only the first cc| superdiagonals different 
from zero, therefore (1) of this Proposition leads to the stated result. 

(3) It follows from Yo = I. 

(4) It is a direct consequence of Proposition |2.21|(3) 

(5) From the string equation ( |2.13[ > and the persymmetry property ( |2.9| > we get 


Oil 


qY«G =r|GY a 
= G T qY a , 

and using the Gauss-Borel factorization ( |2.5[ ) we obtain 

TiY a S -1 H(§ _1 ) t =(S -1 H(§^ 1 ) t ) T TiY 
= (§) _1 H T (S T ) _1 qY 

so that 

%T a S~ ] H =(S- 1 H(§- 1 ) t ) T qY a 

=H T (S T )- 1 T 1 Y a §t 
=H T (S T )~ 1 YT n §t ) 

consequently C a H = H T C^ and C a = I l C a H 1 follows. In particular, CoH = 11 C Q ' 

C 0 = = H T (S T r 1 q$ t H^ 1 , 
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and, as we have a lower triangular matrix on the LHS and a upper triangular matrix on the RHS, 
the only option for Co is to be a diagonal matrix and the second relation Co = r| is proven. Finally, 
notice that 

C* 1 = 

= sr q V 1 ^ 1 

and using Propositions |4.12] and |2.2l] we conclude the third relation. 

( 6 ) Observe that 

C a =(H T )- 1 §pS^ST a S ~ 1 


C.6. Proof of Proposition |3.13[ 

Proof. <= We need to prove that £(z“ L) = | ex'| + m for all a' G Z D ; i.e., that the transformed Newton 
polytope {ex'}+ NP(L) intersects with the faces of Conv([|a'| + ml ) for each a' <E Z D . From Proposition |2.16| 
we know that there exists an orthant labelled by cr G 2 Z ° such that a 1 G (Z D ) a . By hypothesis, for each 
a G 2 Z ° we can find G (M d )cf n NP(L) so that |cxo-l = nr and we can ensure that |a' + a a | = |a'| + |a CT |, 
Vcd G (Z D ) (J . Finally, for each orthant (M D ) a we have such multi-index, consequently in the transformed 
Newton polytope {a 7 } + NP(L) we have a non empty intersection with the faces of Conv([m + cx'|]) for all 
a' G Z. _ 

=> Suppose that there is an orthant (M D ) (T in where we have no points in the Newton polynomial of 
longitude m then for any a' G (Z D ) ff in that orthant we have that £(z“ L) < la'I + £(L). □ 


C.7. Proof of Proposition |3.19[ 

Proof. => Let Li and L 2 be two nice Laurent polynomials. Then, for each subset cr G 2 Z ° there are multi¬ 
indices (cxOa G NIL Li) n (M D ) a of longitude = f f Li), i g {1,2}. Hence, recalling ( 3.11 [ l we see that 

+ (ai)ff G NP(LiL 2 ) and moreover, as both belong to the same orthant, |(ai) ff + (c^o-l = l(cxi)o-l + 
l(a 2 )(j| = f(Li) + £(L 2 ). Therefore, (a^o-T G NP(L 1 L 2 ) nConvftffLi) + £{L 2 )]) Cl (M D )o- f 0/ and 


recalling Proposition [3T4] the result is proven. 

<= Let Li and L 2 be two Laurent polynomials being at least one of them not nice, assume that Li is not 
nice. This means that there exists at least an orthant where all the multi-indices in NILL]) n (M D ) (y fulfill 
| (cci) o-| < £ (Li). Thus , when we consider the product Li L 2 and we look in this orthant (M D ), y we will have 
multi-indexes of the form (ai) CT + (a 2 )o- with |(ai) a + (a 2 )al = l(ai) CT | + I (CX 2 ) cr I < £(Li) + £(L 2 ), and the 
equality is never saturated. The sum of multi-indexes of different orthants will not help in achieving the 
longitude £(Li) + £(L 2 ). □ 


C. 8 . Proof of Lemma lL2l 


Proof. We will show that the | [k + 1 ] | columns in this matrix are actually linearly independent. For that to 
be true it is enough to show that the only solutions for the equation (x[i] (T)) ^ k ] X = 0, with 


X = 


/ X a‘ ,c+1) \ 


\ x a (k+l) / 


G C |[k+1]| , 


is X = 0. First, let us consider the equation (T a )f lc ]p c+ i]X 
the previous equation will hold as long as every X ( k +i) 

a i 

holds for (Tq 1 )[k][k+l]X = 0, it implies that every X^k+i) 


= 0 from the Definition 

ik+1) 


2.20 


we can deduce that 


such that a- a > 0 is equal to zero. The same 


with ct| k+1) < 0 must equal zero. Therefore, for 
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both equations to hold, we need that every X^, k 11 ; having a- k Q 1 f 0 must vanish. In other words, only 
those X^fk+i] having = 0 can be different from zero. But our request is for all a G {1,2,..., D} 

which implies that the only nonzero component of X is the one for a = (0,..., 0) T 0 [k + 1] and therefore 
X = 0. □ 

C.9. Proof of Lemma [431 


Proof. For the first line just consider the differentiation of the deformed LU factorization (4.41 in the 
form 

(C.l) 

which gives 


-i\t 


S(t)Wo(t)G = H(t)((S(t))~ 1 ) 


9S S^ 1 + SY a S _1 = ( - H (V 1 ) 


t 9$+_! 


3t a V9t a 

Split it into upper and lower triangular to get 

3S_i _i /_i\ / 3H 


3t a 


FF 


^-S _1 +SY a S _1 = (sr a S~ ] ) = H _1 — H f§ _1 V H -1 ), =-fsY a S 

3t a v )> \3t a V / 3t a / 3t a V 

While for the second one consider the differentiation g^— of the deformed LU factorization ( |4.4| in the 
form 

(C.2) 

which gives 


3S 




$(t) (W 0 (t)) t Gt = (H(t)) t ((S(t))- 1 ) t . 

t 3S+ 


i i / 3FF 

-S~ ] + Sr a S~ ] = (-= 

3t_„ V 3t 




Now, using the same splitting technique as before 


3S 

3t o 


§-! + §ya§-1 = ( $ya$- 




/s 


- — I SY a S 

3t_ a V /< 




□ 


C.10. Proof of Proposition |4.7[ 

Proof. (1) By differentiation in Definition [43] we get 


K*r s_1 + sr “ s_1 ) Wl ’ 

otcc V 9t a / 


3W 


3S 


§ _1 + §Y«§- 1 jW 1) 

tH-’j’-Mt ( S -) t ,W 2 . 




3t-a v 3t_ a 

3W 2 / 3Ht 


3t_ a 


(2) From Definition 


2.22 


Consequently, fr om L emma 4.3 and its proof we get the result. 

we know that the perturbed Jacobi matrices reads: J a (t) := S(t)Y a (S(t)) 1 


and Ja(t) := §(t)Y a (§(t)) \ therefore by differentiation we get 


3Jc 

3t„ 


95 S" 1 1 

° ) ) OC 


3t 


3 Jc 


3t_ c 


r 3S 

3t-« 


0 }) a 


and with Lemma 1431 we obtain the stated result. 
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(3) As a compatibilitiy condition we can write 


3 2 W 1 

3t«9t a / 



3 2 W, 

3t«'3t« 


and therefore we deduce the zero-curvature equation. For the rest of the expressions we proceed 
similarly. 


□ 


C.ll. Proof of Proposition |4.16[ 

Proof. In order to prove the first equation we notice that 

Sfa^W! =(3 Q ,bS + ST Q T b )W 0 , 

SaSbW! = (3 a 3 b S + 3 a ST b + 3 b ST a + ST a T b )W 0 


and therefore (3( Q)b ) — 3 a 3 b )Wi = — (3 a S(T b )> + 3 b S(T Q )>)W 0 + 1W 0 so that 

(3(a,b) - 3a3b + V Q)b + V b ,a) Wi G (W 0 . 

On the other hand, it is obvious that 

^3(a,b) — 3 Q 3 b + V Qjb + V b)Q )w 2 G u. 

Now, we apply Proposition |4.13| with 

Ri = ^9(a,b) — 3 a 3 b + V Qib + V bjQ j Wi, i = 1,2, 

to get the first result. In order to prove the second result proceed in the exact same way s tarting with 
Wi, deriving with respect to 3 j_ Q! bj and 3_ Q and finally use the results in Proposition 


(tW 0 )t. 


4.13 


that involve 

□ 


C.12. Proof of Proposition |4A8j 

Proof. We begin with 

(3(a,b,c) - SaSbScJtW,) =(3 Q)b , c S - 3a3 b 3 c S - 3 a 3 b ST c - 3 b 3 c ST a — 3 c 3 a ST b 

- 3 a ST b r c - 3 b ST c T a - 3 c ST Q T b )W 0 , 

and take into account S = I + |3 (1 * + |3 (2 * + • • •, being (3'■ k ' the k-th subdiagonal of S, we can write 

(9(a,b,c) - 3a3b3c)(Wi) = ( - 3 a 3bP ci:1 (T c )> - 3 b 3c(3 m (T a ) > - 3 c 3 a (3 (1) (T b )> 

- (1 - 6o,b+c)3a(3 (2) (T b )>(T c ) > - (1 - 6 0) a + c)3b(3 (2) (T c )>(T a ) > - (1 - 6o,a+b)3c(3 t2) (T Q )>(T b ) > 

- (1 - 6o ) b +c )3a(3 (1) (T b )>(T c ) > - (1 - 6o,a + c)3 b (3 (1) (T c )>(T a ) > - (1 - 6o,a + b)3c(3 m (T a )>(T b ) > )Wo + (W 0 

Now recall 


3aW! = ((T q )> + (3(T a)>)Wo + tw 0 

and use it in every term of the third line of the previous expression to obtain 

(3(a,b,c)-3a9 b 3 c + (l -6 0 ,b+c)3a(3 (1) (T b ) > 3 c + (l -5 0 ,a + c)3b(3 (1) (T c )>3 Q + (l -6o,a+b)3 c (3 (,) (T a ) > 3 b )(W 1 ) = 
(-3 a 3 b (3 t1] (T c ) > -3 b a c p (1 )(T a ) > -3 c 3 a p (1) (T b ) > -(l -6o,b+c)3a(3 t2) (T b )>(T c ) > 
-(l-6o,a+c)3 b |3 t2) (T c ) > (T a ) > -(l -6o ; a + b)3 c (3 (2) (Ta) > (T b ) > + (l -6 0 . b+c )3a(3 (1) (T b ) > (3 ll) (T c ) > 

+(1 - 5o,a+c)3 b |3 (1) (T c )>|3 (1) (T a ) > + (1 - 6o,a+b)3 c p (1) (T a ) > P (1) (T b ) > )W 0 + OVt . 

From where the result follows. □ 
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